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A general theory of analytic modulation systems is developed where 
the transmitted signal is of the form a(t) = Re \e lWct f(z(t))\. Here f(z) 
is an analytic function (modulation law), and z(t) = x(t) + iy(t) is the 
analytic baseband signal whose real part x(t) is a bounded bandlimited 
signal of spectral support [— ft, ft] which is assumed to have a bounded 
Hilbert transform y(t). It is shown for a large class of \z(t)\ and modu- 
lation laws that z(t) may be recovered using a receiver incorporating 
the inverse function of fas a detector with appropriate pre- and post- 
detection filtering. The theory also shows that in the procedure for 
factoring certain positive bandlimited signals, an approximate Hilbert 
transform operator (bandlimited) may be used. A related result is that 
signals subjected to logarithmic companding (one-sided) and filtering 
may be recovered by a non-iterative method. 

I. INTRODUCTION 

In 1962, Bedrosian proposed a modulation system called single side- 
band phase (or frequency) modulation.* (See Ref. 5.) The modulated 
signal is of the form 

Re { e <"<a^+*u)+iyU))| = e -y(t) cos ( Wc£ + x ( t )) 

where x(t) is the "baseband" signal, y(t) is the Hilbert transform of x(t), 
and u c is the carrier frequency. The special relation between the am- 
plitude modulation e~ y{t) and the phase modulation x(t) results in the 
modulated signal having no spectrum in the interval (— a) c ,co c ); i.e., the 
amplitude modulation removes the lower sideband. However, the 
spectrum is still infinite in extent. We adopt the terminology "single 
sideband exponential modulation" (SSBEM) for this system. 

Bedrosian pointed out that SSBEM was compatible with conventional 
FM receivers, and suggested that the single-sideband system might offer 



* K. H. Powers received a U. S. Patent (No. 3,054,073) on such a system shortly before the 
appearance of Bedrosian 's paper. See Voelcker. 21 
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some savings in transmission bandwidth over the conventional system. 
However, since filtering operations could radically alter the zero crossings 
of the modulated signal, it was not clear to what degree one could 
maintain compatibility and at the same time realize some saving in 
bandwidth. • 

Others 2,10 - 17 have compared the spectral distribution of single side- 
band and conventional frequency modulated signals for the cases of si- 
nusoidal and Gaussian noise modulation. They have shown that the 
"effective" bandwidths, as measured by central second moments, of the 
single-sideband signals may be greater or less than that of the conven- 
tional FM signal depending on the nature of the modulation. At any rate, 
it is not clear how one would translate these results into relative band- 
width requirements of the two systems, each employing a conventional 
FM receiver. 

Aside from the compatibility question, Barnard 4 has shown that the 
transmission bandwidth requirements of single-sideband exponential 
modulation are, in a strict sense, minimal. He showed that if the mod- 
ulation x(t) belonged to a certain subclass of bandlimited signals with 
spectral support [— ft,ft], that the modulation could be recovered, within 
an additive constant, from a knowledge of the spectral distribution of 
the single-sideband signal in the interval [u; c ,u> c + 12 + e], provided 6 > 
0. This was proved by demonstrating the convergence of an iterative 
recovery scheme. 

Here we consider a class of single-sideband systems wherein the 
modulated signal is of the form 

s(t)=Re\f(z(t))e i °* t ] 

where f(z), the "modulation law", is an analytic function and z(t) =x(t) 
+ iy(t) is the "analytic signal" of which, say, the real partx(t) is the in- 
formation to be transmitted. We suppose that x(t) is bounded and 
bandlimited with spectral support [-ft, ft] and that s(t) is transmitted 
over a channel whose transmission function is the Fourier transform of 
an absolutely integrable function (impulse response) and is equal to unity 
over (a> c ,a> c + a). It is shown under fairly weak conditions on f(z) and 
z(t) that x(t) may be recovered (by a relatively simple non-iterative 
method) from the received signal, provided a > ft. 

The gist of the method can be grasped by considering periodic signals; 
e.g., 

~ " (X = 0, say) 

z(t) = tx k e ikt . 

i 
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We suppose that sup | z (t ) \ = m, and f(z ), the modulation law, is an- 
alytic for |2| < m,(/(0) = 0, say) 

f(z) = Y.a k z k , \z\ <m. 

1 

Then setting w(t) = f\z(t)\ we have 

w(t) = £ W^e'*' 
1 

where W k depends only on those Xj and a/ for which l<j <k. This sort 
of dependence allows us to determine x(t) from a bandlimited version 
of w(t), say, 

u>n(t) = L W k e ikt , 
1 

by what amounts to reversion of power series, provided /'(0) ^ 0. We 
have 

z = <f>(w) = £ bftitf*, for |u;| sufficiently small. 
1 

(61 = ( fll )-i = {/'(0)}- 1 ) 

Assuming that the series converges when w is replaced by w„(£) we 
set 

Zn(t) = <t>\w n (t)) = T.b k \w n {t)\ k 

1 

and then by formal composition of the power series find that 

Znit)=*t.X k e ikt + £ c*e'*'. 

1 n+i 

So the first n Fourier coefficients of z(t) and z n (t) agree; i.e., under the 
stated assumptions, x(t) can be recovered by bandlimiting <{>\w n (t)\, 
where is the inverse function, z = <l>(w), and w n (t) is the partial sum 
of the Fourier series of w (t ). 

The simplicity of this procedure owes to the fact that z (t) has a one- 
sided spectrum and the analytic modulation law f(z) then gives a 
function w(t) = f\z{t)\ which also has a one-sided spectrum. Since z{t) 
contains no negative-frequency components, the usual difference terms 
do not appear; i.e., the spectrum of w(t) in the frequency interval [0,cv] 
depends only on the spectrum of z(t) in the same interval. 

The recovery procedure is not so transparent for more general band- 
limited signals x(t). First of all, filtering w(t) with a filter whose trans- 
mission function is unity over [0,a] and zero for frequencies greater than 
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(3 (i0 > a) will give a function w ai p(t), analogous to w n (t), which may 
differ considerably from w(t). It may be that <f>\w ai n(t)\ does not have 
a one-sided spectrum; i.e., <t>\w ail j(T)\, t = t + iu, is not analytic in the 
upper half-plane u > 0. Indeed w ailj {t), (— « < i < °°), may not even be 
in the domain of definition of </>; i.e., </> may have a natural boundary 
beyond which it cannot be extended. Even if (f>\w a ^(t)\ does have a 
one-sided spectrum, the function <t> need not have a power series repre- 
sentation over the range of w a ^{t) so that one cannot use convolution 
arguments to show that the Fourier transforms of (f>[w(t)\ and (f>\w ai p(t)\ 
agree over [0,a\. This particular problem is met by using generalizations 
of the Paley- Wiener theorem. 

The problem arising when <j>\w„ ,#(£)} does not have a one-sided spec- 
trum (w a> p(t) not in the range of the inverse function) is met by imposing 
restrictions on z(t), namely that for sufficiently large u, the range of z(t 
+ iu) is sufficiently small that w(t + iu) will be in the range of the inverse 
function. This implies that w a ,^{t) may be filtered (with a Poisson filter) 
to obtain w ai/ i(t + ib), which for sufficiently large b will be in the range 
of the inverse function. One can then obtain z (t + ib) and then use in- 
verse Poisson filtering to recover z(t). 

Although one could conceivably recover z(t) from w„ t p(t) by other 
procedures when ^{w^uir)] is not analytic for u > b, the method here 
avoids any decision process and gives a simple receiver model incorpo- 
rating the inverse function and (possibly) a Poisson filter with its 
(bandlimited) inverse and appropriate low-pass filter. 

Generally speaking, given f(z) and the channel transmission function, 
one can design a receiver which will work for a certain subclass \z(t)\ of 
signals. Or given f(z) and a fixed receiver design, one may ask for the 
minimum bandwidth channel required for transmitting a given subclass 
of signals. In this connection, some estimates are given for the bandwidth 
requirements of "compatible" single-sideband exponential modulation 
with \z (01 all functions of spectral support [0,1] satisfying sup \z{t)\ < 
m. In a recent work, Werner 23 has considered the same problem for z(t) 
in Z/2 and gives upperbounds in terms of the Z/2-norm of 2. 

There are rather dramatic mathematical simplifications in the de- 
tection theory when the signals \x(t)\ are restricted to be of the band-pass 
type, allowing radical changes in the system design. 

The theory also shows that the factorization of certain positive band- 
limited signals can be effected with an approximate Hilbert transform 
operator acting on the logarithm of the signal. A related result pertains 
to the signal recovery problem considered by Landau and Miranker 11,12 ; 
viz., there is one companding function ("log") for which the signal can 
be recovered by a non-iterative method. 

Another interesting consequence of the general theory is the fact that 
for n arbitrary numbers a*., k = 1,2, . . . , n, there exists an integer v > 
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n and corresponding numbers ctfe, k — n + 1, . . . , i>, such that the poly- 
nomial 

P r (z) = l + ta k z k 
i 

is zero-free for \z\ < 1. 

Although the general theory is interesting from a mathematical 
viewpoint, it would appear that the practical interest in analytic mod- 
ulation systems, other than the linear system, is limited to SSBEM, i.e., 
to the case f(z) = e z (or e lz ). In this case one can trade bandwidth for 
simplicity of detection. However, the trade-off is attractive only for 
moderate amplitudes of z(t) where SSBEM offers an interesting alter- 
native to other systems employing envelope detection. 

It should be noted that the method here is naturally confined to 
bounded bandlimited signals z (t), since otherwise we would require both 
f(z) and its inverse 0(u>) to be entire functions, i.e., w = f(z) = a + bz. 
The exception would be the band-pass case where / could be an entire 
function and <i> replaced by an equivalent polynomial (see Section 5). 

Of course, the theory here has to be extrapolated to practice with the 
appropriate "epsilons"; i.e., an analytic modulation law/(z) can only be 
approximated within t\ over the disk \z\ < m, and an analytic signal z(t) 
having one-sided spectrum can be realized in practice within 62, and the 
impulse response of an ideal filter can be approximated (in Li) within 
€3, etc. Then the continuity of the overall transformation may be used 
to bound the errors. 

In order to deal rigorously with "communication type" signals which 
do not have ordinary Fourier transforms a considerable amount of pre- 
liminary mathematics is required. However, one can follow the theory 
assuming that the signals are either periodic or have ordinary Fourier 
transforms, with one cautionary note in mind. Abrupt bandlimiting 
operations (spectral projections) such as convolution with sin t/irt are 
not permissible (not defined) for the general signals of interest. 



II. PRELIMINARIES 

A measurable function g{t) is said to belong to L p (— °°,«>) abbreviated 
hereafter as L p , (1 < p < °°) if 



X- lg{t)]i 



\Pdt < 00. 
The L p - norm of g is defined by 

llsllp = U_~ \g{t)\Pdt\ 
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and if a is a scalar 

Hag lip - |«| • Hi? lip- 

If \g{t) | is uniformly bounded, with the possible exclusion of a set of 
measure zero, g is said to belong to L m and the norm of g is 

||^ || co = essup \g(t)\ 

t 

where "essup" over t is the essential supremum of |g(£)|, which is the 
infimum of numbers M such that 

\g(t)\ < M for almost all t. 

We will be mainly concerned with continuous bounded functions g(t) 
in which case 

||*||. «= sup \g(t)\. 
t 

For 1 < p < a>, the L p norm satisfies the triangle inequality 

Wgl+gA^ 11*1 lip + Ifslp (i) 

which for any sequence of numbers ak satisfying 2 1 a^ | < oo and se- 
quences of functions gk such that ||** || < M, leads to 

WLakgkWp * ZWklhkWp- (2) 

There are functions which belong to L p for only one value of p. 
However, it is easy to see by considering the set where \g(t)\ < 1 and the 
set where \g(t) \ > 1, that if* belongs to L r and L s where 1 < r < s, then 
g belongs to L p for every p satisfying r < p < s. For example, the func- 
tion sin tit belongs to L p for every p > 1. 

Associated with the space L p is the conjugate or complementary space 
L q where 

i + A-L 

P Q 

For functions in complementary spaces we have Holder's inequality 



S- 



g(t)h(t)dt 



< WgWJhL, l<p<oo, (3) 



which for the case p = q = 2 is the familiar Schwarz's inequality. 

In connection with Holder's inequality we note that the norm of a 
function may be equivalently defined as 



||*|| p = sup f h{t)g{t)dt 



\\h\\ q = l,q = p/(p-l). 
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A convolution kernel K in L\ carries L p into L p . We have 

K9g(t)= f" g(s)K(t-s)ds i mL r p 
J-<* KmLi 

and 

\\K ® g\\ p < \\KU\gWp (4) 

which amounts to a generalization of (2) to weighted sums of translates 
of g. The convolution integral is not in general defined for each t unless 
K belongs to the conjugate space of g. In general the convolution is de- 
fined as the limit in L p of 



gm(t) = §~j{s)K m {t- 



s)ds 



where K m is a sequence of bounded functions of Li (hence K m belongs 
to L q ) satisfying 

lim 



im f " \K(t) - K m (t)\dt = 0. 

— »OD %J — CD 



2. 1 The Fourier transform on L p 



A function g in L p has an ordinary Fourier transform, provided 1 < 
p < 2, (a theorem of M. Riesz, cf. Ref. 20) in the sense that 



grM = C_ T g(t)e-^dt 



converges in norm as T —*■ » to a function g(u) belonging to the com- 
plementary space L q , i.e., there exist a function £ in L q such that 

lim \\g-g T \\ q = 0. 

However the Fourier transform on L p , 1 < p < 2, does not carry L p into 
all of L q except in the case p = 2. In particular, the Fourier transform 
of a function g of L\ is a continuous function. Furthermore, (the Rie- 
mann-Lebesgue Lemma) 

lim f g(t)e~ iwt dt = 

a, —±oo c/-« 

for g in L\. Unfortunately, there is no simple description of functions 
g(co) which are the Fourier transforms of functions g(t) of L\. A useful 
sufficient cordition is that^(co) belong to Li and have a "derivative in 
L 2 " (meaning only that£M is the integral of a function of L 2 denoted 
by dg/dw). The sufficiency of this condition may be seen by writing 
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C"\g(t)\dt = f " , }. r \{t + ia)g(t)\dt (a>0) 

and then applying Schwarz's inequality and Parseval's theorem. The 
result is (choosing the best value of a ) 



is: 



\g(t)\dt 



< Wh- 



(5) 



In obtaining this result we use the fact that § in L 2 and dg/du in L 2 imply 
that£(a>) tends to zero at ±». (By Schwarz's inequality, gdg/du belongs 
toLi, so |£(<o)} 2 is absolutely continuous and tends to limits at ±«. The 
limits must be zero in order for £ to belong to L 2 .) 

2.2 Bounded functions whose Fourier transforms vanish over certain sets 

It is not necessary to attempt to define the Fourier transform of a 
bounded function g(t) in order to give precise meaning to the statement 
that the Fourier transform ofg(t) vanishes over some open set E. This 
can be done in a way which is consistent with the ordinary Fourier 
transform, should it exist, of g(t) vanishing over E. Here we restrict E 
to be the union of a finite number of disjoint open intervals. 

Definition: The Fourier transform of a bounded function g is said to 
vanish over E if and only if 

(i) CJj(t)h(t)dt - 

for all hinLi whose Fourier transforms satisfy 

{ii) r\(w) = C~ h(t)e~ i * t dt = 0, u<£E. 

This definition has its logical basis in Parseval's formula for functions 
of L 2 . The bar over h in (i) denotes the complex conjugate of h. It is 
readily verified that (i) may be replaced by 



(Hi) f" g(t)h(-t)dt = 



which is more directly applicable to convolutions. That is, if the Fourier 
transform of g vanishes over E we have 



ig ® h)(t) = f " g(s)h(t - s)ds = 



(6) 



for all h'mLi whose Fourier transforms vanish outside E. 

We also note that in case the set E is symmetric with respect to the 
origin, h{t) in (i) may be replaced by h(t). [See Ex. 4 below.] 
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We say that the Fourier transforms of two bounded functions g\ and 
g 2 agree over E if and only if the Fourier transform of {gi — g 2 ) vanishes 
over E. We also say that g has spectral support E c (a closed set), or the 
spectrum of g is confined to E c , meaning that the Fourier transform of 
g vanishes over the complement of E c . 

The following are some elementary consequences of the definition. 
The proofs are left as simple exercises. It is understood throughout that 
g, g\, and g 2 are bounded functions. 

Example 1. Suppose all the intervals composing E are finite and K{t) 
is any function of L\ whose Fourier transform K(to) satisfies 

#(«) = 1 for wGfi. (7) 

Let 

g 2 (t) = j'^g l (s)K(t-s)ds. (8) 

Then the Fourier transforms of gi and g 2 agree over E. 

Example 2. If the Fourier transform of g(t) vanishes over («,0), then 
the Fourier transform of e lXt g{t) vanishes over (a + X,/3 + X). 

Example 3. If the Fourier transform of g\ vanishes over E\ and the 
Fourier transform of g 2 vanishes over E 2 , then the Fourier transform 
°f (Si + 82) vanishes over E\ C\ E 2 . 

Example 4. If the Fourier transform of g vanishes over E, then the 
Fourier transform of the complex conjugate g vanishes over E^~\ where 
E^~) denotes the reflection of E with respect to the origin. 

Example 5. If the Fourier transform of g vanishes over E, then the 
Fourier transform of Re \g\ (or I m \g\) vanishes over E C\ E { ~\ 

Note: E (~) E ( ~^ may be the null set. However, if E is symmetric with 
respect to the origin, E = E { ~\ Hence a class of functions whose Fourier 
transforms vanish over a set E which is symmetric with respect to the 
origin is essentially a class of real-valued functions, since the real and 
imaginary parts of the functions separately belong to the class. 

Example 6. (Reproducing Kernels) Suppose the spectrum of g is con- 
fined to a set E c consisting of n finite disjoint closed intervals. Let K(t) 
be any function of L\ whose Fourier transform K(co) satisfies 

K(w) = l, u>EE c . (9) 

Then for almost all t we have 

g(t) = £"g(s)KU-s)ds. (10) 
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(Set g\ = K ® g and show that (g\ — g) is orthogonal to all of L\.) 

Note: The qualification "almost all t" arises because the definition of 
a function g on a set of measure zero is irrelevant to the condition for its 
Fourier transform to vanish outside E c . However, (K ® g)(t) is a con- 
tinuous function of t and s'o we will adopt the convention that a function 
such as g in Ex. 6 is continuous. 

We note further that the condition that K in Ex. 6 belong to L\ can 
be relaxed in case g belongs to L p for some p satisfying 1 < p < °°. In this 
case we can take K(a>) = for a> (£ E c . It is sufficient to prove this when 
E c is a single interval [— fi,12] and this has been done (Ref. 14). 

2.3 The Paley-Wiener Theorems for Loo 

There is an important connection between functions whose Fourier 
transforms vanish over a half-line and functions analytic and of expo- 
nential type in a half-plane. The following theorems are extensions to 
La, of the classical "one-sided" and "two-sided" Paley-Wiener Theo- 
rems 18 for Li. 

Theorem 1. The Fourier transform of a bounded function g vanishes 
over (— »,a) if and only ifg(t) is the boundary value of a function g{r), 
r = t + iu, analytic in the upper half-plane u> and satisfying 

sup \g(t + iu)\ <e~ au sup \g(t)\ for u > 0. (11) 

t t 

There is the analogous theorem connecting functions g(t) whose Fourier 
transforms vanish over (/?,<») and functions g(r) analytic in the lower half 
plane. The specialization of Theorem 1 to functions whose Fourier 
transforms vanish over (-«>,«) and (/3,°°) is the following. (We assume 
that — <» <a<j8<°° and according to the convention above qualify g 
to be continuous.) 

Theorem 2. The Fourier transform of a continuous bounded function 
g vanishes outside [«,/?] if and only ifg(t) is the restriction to the real 
line of an entire function g(r), r = t + iu, satisfying 

sup \g(t + iu)\< e~""sup \g(t)\, u>0 
' ' (12) 

<e~0"sup \g(t)\, u<0. 

t 

These theorems are essential to the theory of single-sideband systems 
for bounded signals which do not have ordinary Fourier transforms. 

Actually we do not need a uniform bound on the rate of growth (decay) 
of g{t + iu) to infer that the Fourier transform of g(t) vanishes over 
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(—w, a). In fact, an asymptotic bound implies a uniform bound. 
Theorem 3. If g(r) is analytic in the upper half-plane and satisfies 

sup \g(t + iu)\ < co for u > 

t 

then the asymptotic estimate 

sup \g(t + iu)\ = 0\e~ au \ asu — *> 

implies 

sup \g(t + iu)\ < e~" u sup \g(t)\ for u > 0. 

; t 

Proofs of Theorems 1, 2, and 3 are given in Appendix A. 

We note the following corollaries of Theorems 1 and 2, concerning the 
Fourier transforms of products. 

Corollary 1. If the Fourier transform of gi uanishes over (—«>,«) and 
the Fourier transform of g-2 uanishes over (— <x ',l3), then the Fourier 
transform of g\g2 vanishes over (— °°,« + fi). 

Corollary 2. If the Fourier transform of g\ vanishes outside [«i,/Si] and 
the Fourier transform of g% vanishes outside [a^foL then the Fourier 
transform of g\gi vanishes outside [«i + «2,/3i + fo]- 

2.4 Terminology 

Functions whose Fourier transforms vanish outside a finite interval 
are called bandlimited functions. Generally, we think of the interval 
centered at the origin and refer to bandlimited functions also as low-pass 
functions. 

Functions whose Fourier transforms vanish over an interval centered 
at the origin are called high- pass functions, and functions which are both 
high-pass and low-pass are called band-pass functions. 

Functions (signals) whose Fourier transforms vanish over a half-line, 
usually (— °°,0), are generally called analytic signals. 

2.5 The Hilbert transform and the analytic signal 

We would like to map the space of real-valued bounded signals x(t) 
of spectral support [ — S2,i2] into the space of complex-valued bounded 
signals z(t) of spectral support [0,8]. We would like the mapping to be 
linear and also have the property that translates of x map into translates 
of 2, so that no "time stretching" is involved. The usual way of doing this 
is to take 

z(t)=x(t) + iy(t) (13) 

where y = x, the Hilbert transform of x. 
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The Hilbert transform is defined by 

i(t)=--f"— ds (14) 

7T J-<*> t — S 

where the cut in the integral sign indicates a Cauchy principal value at 
s = t. The difficulty we encounter is that an arbitrary bounded band- 
limited function does not have a Hilbert transform so we have to restrict 
x somehow. 

We may regard the Hilbert transform as the limit as a -*■ of convo- 
lution transforms 



l a (t)= C™ x(s)K a (t - s)ds 



(15) 



with regular kernels K a given by 



/t a (t) = --i~-£,a>0. (16) 

Now K a (t) belongs to L p for every p > 1 and has a Fourier transform 
R a (a>) given by 

£ a (w) = -i(Bgn«)e-«M. (17) 

Thus if x(£) has a Fourier transform x(u>), the Hilbert transform x{t) has 
a Fourier transform given by 



S- 



x{t)e- iuJt dt = -i(sgn «)£(«) (18) 



and consequently the Fourier transform of z(t) as defined in (13) van- 
ishes over (—<»,0). Now 2(0 is the boundary value of the function z(t), 
t = t + iu, defined by 

*( r )«i C"^-ds, u>0. (19) 

IT •/-» T — S 

In case x does not have an ordinary Fourier transform, but has a bounded 
Hilbert transform x, the function z(t) is bounded and analytic in the 
upper half-plane and according to our definition and Theorem 1, the 
Fourier transform of z(t) vanishes over (— <»,0). Also if the Fourier 
transform of x(t) vanishes outside [—0,12] then the Fourier transform 
of z(t) vanishes outside [0,Q]. 

The subclass of bounded functions which have bounded Hilbert 
transforms does not have a simple alternate description. In order for x(t) 
to have a bounded Hilbert transform it is sufficient that x(t) have a 
bounded derivative and a bounded integral (Ref. 15). If the Fourier 
transform of x vanishes outside [— fl,ii] then x (t) has a bounded deriv- 
ative ("Bernstein's Theorem," cf. Theorem 11.12, Ref. 6) satisfying 
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sup|x'(*)| <flsup|x(t)| f (20) 

t I 

and there are various ways to restict x(t) to have a bounded integral. For 
example, we may assume that g{t) is an arbitrary bandlimited signal and 
set 

x(t)=g'(t). (21) 

(Also any high-pass function has a bounded integral (Ref. 14).) Then we 
can determine g(t) within an additive constant from the real part of z(t). 
Assuming that the Fourier transform of g vanishes outside [-ft, ft] we 
have the inequality (Theorem 11.4.3, Ref. 6) implied by (21), 

sup|x(0| <flsup|g(0|. (22) 

t t 

An interesting subclass of bandlimited functions which have bounded 
Hilbert transforms are the band-pass functions. For these functions, 
there are equivalent Hilbert transform kernels which belong to L\. If the 
Fourier transform of x(t) vanishes outside the intervals [— ft,— rft] and 
[rft,ft] where < r < 1, then x has a bounded Hilbert transform satisfying 
(Ref. 7) 



sup |x(t)| ^ 



2, 1 

A + - log - 

it r 



sup|x(0| (23) 



where A < A/it and 2Ar cannot be replaced by a smaller number (i.e., as 
r approaches zero). 

Other subclasses worthy of mentioning may be generated by convo- 
lution transforms on L^. Thus if |g| < M and k is a kernel in L\ which 
has a Hilbert transform also in L\ then the class of functions of the 
form 

x(t)= f" g(s)k(t - s)ds (24) 

have bounded Hilbert transforms given by 

x(t)= C~ g(s)k(t-s)ds. (25) 

For example, x{t ) may be the output of some crude sort of band-pass 
filter (like an a-c amplifier) as would be the case for 

k{t) = ae~ at -be~ bt , t>0 

= 0, t < (26) 

where a > b > 0. Then 

r-arc.-Hft--^-;-* <°- fc) ';. ,27) 

Jo a + ico + 100 (a + iu)){b + Ha) 
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5- 



and 

me-'»>dt = <",- ;>H. (28) 

(a + tu;)(o + ico) 

It follows from (28) and (5) that k belongs to L\. 

Another sufficient condition for a bandlimited function x(t) to have 
a bounded Hilbert transform is the requirement that x belong to L p 
where 1 < p < «. Then if the Fourier transform of x vanishes outside 
[-fl.fl] we have (Ref. 14) 

s in Sl(t — s ) 
ir(t -s) 
Then x has a Hilbert transform given by 

1 — cos S2(£ — s) 
Tr(t -s) 



C" , - sin 9.{t — s) , 
x(t)= ( x(s) -f ^-ds. (29) 

J-a, Tr(t — S) 



m- r"»(.) 1 "" Bff r* , A (30 

*/-<= 7r(t — s) 

and thus by Holder's inequality 



w>i*i*Ip{£ 



1 — cos Qt 



TTt 



"dt 



\hl 



where - + - = 1. (31) 
P Q 

However, the condition that x belong to L p (1 < p < <*>) is not a sat- 
isfactory condition for communication signals. 

Hereafter we will suppose that x{t) is so restricted that z (t) = x(t) + 
iy(t) is a bounded function whose Fourier transform vanishes outside 
[0,12]. We should note that this assumption does not imply thaty is the 
Hilbert transform of x (even within an additive constant). That is, the 
assumption does not imply that the integral 



X 



T x(t)-X(0) 

at 

-T t 



tends to a limit as T — *■ °°, so in effect we are allowing some functions that 
are not of the form z{t) = x(t) + ix(t). In case we further restrict z(t) 
to be a bounded function whose Fourier transform vanishes outside 
[rfi,ft] where < r < 1, then we can assert that z (t) is of the form x(t) + 
ix(t) where the Fourier transforms of x and x vanish outside [— ft,— rfi] 
and [rft,fi]. Conversely if x is a real-valued (band-pass) function whose 
Fourier transform vanishes outside these two intervals, then x(t) has 
a bounded Hilbert transform x(t) and \x(t) + ix(t)} is a bounded function 
whose Fourier transform vanishes outside [rQ,Q]. In other words, we can 
always regard any real-valued band-pass signal x(t) as the real part of 
an analytic signal z(t). This is a special case of a representation theorem 
for high-pass signals (Ref. 13). 
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III. MODULATION AND EQUIVALENT BASEBAND TRANSMISSION 

We suppose that z(t) is a bounded bandlimited signal with spectrum 
confined to [0,12]. Now let f(z) be a function which is analytic over a re- 
gion which includes the disk \z\ < m, where m = sup \z(t) | . We take f(z) 
as a modulation law and generate 

w(t) = fa{t)\ (32) 

which is the boundary value of a function w(t) bounded and analytic in 
the uhp. Hence the Fourier transform of w(t) vanishes over (— ro ,0). 
Generally, the modulation process also includes translation of the 
spectrum of w(t) leading to a transmitter output 

a(t) = Re \w(t)e^tl (33) 

where u) c > is the carrier frequency. The Fourier transform of a(t) then 
vanishes over (— o> c ,co c ) and hence ait) is called a single-sideband signal, 
although the upper side-band may be infinite in extent. 

In conventional single-sideband amplitude modulation (SSBAM) the 
modulation law is the linear law, f(z) = z, in which case the Fourier 
transform of a(t) vanishes outside the intervals [w c ,oj c + ft] and [— a> c — 
ft, —<0 C ], so that the bandwidth required for transmitting a{t) is (counting 
positive and negative frequencies) 2ft + e where e is an arbitrarily small 
positive number. In other words, a(t) has a reproducing kernel in Li of 
bandwidth slightly larger than 2ft. (Recall that the Fourier transform 
of a function of Ly is continuous.) We will see that the spectral economy 
of SSBAM carries over to more general modulation laws f(z). So we as- 
sume that o(t) is transmitted over a channel (characterized by an Li 
impulse response) which has unity transmission over the frequency 
bands [o> t .,a; c + a] and [— u c — a,— u c ] where a > ft. The transmission 
may be zero outside slightly larger intervals. 

We denote the received signal by a nit) and since its Fourier transform 
vanishes over (— co c ,u; c ) it has a Hilbert transform OR(t). We assume that 
the carrier frequency and phase are known at the receiver so that we can 
form 

w a (t) = e-^loRit) + ia R it)). (34) 

In engineering parlance the real part of w a (t) is obtained by in-phase 
synchronous demodulation of a^(i), while the imaginary part of w a (t) 
is obtained by quadrature synchronous demodulation of <tR(t). The 
Fourier transform of w a (t) vanishes over (— «>,0) and we have 



w, 



i a it)= C w(s)K„{t - s)ds (35) 

where K a (t) is the impulse of an equivalent baseband channel satisfying 
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for some a > 12 

f " K a {t)e~ iut dt = 1, for < a> < a (36) 

\K a {t)\dt<«>. (37) 



X 



Hereafter we will be concerned with recovering z(t), and hence x(t), from 
w a (t) in the equivalent baseband transmission of w(t) as given by 
(35). 

IV. THE INVERSE FUNCTION AS A DETECTOR 

We would like to solve (35) for z(t) where w(t) = f\z(t)\. This is (su- 
perficially) similar to the problem of Landau and Miranker 11,12 where 
w(t) = f\x(t)\ and / is a real function of a real variable, x(t) is a real- 
valued bandlimited function of L^ whose Fourier transform vanishes 
outside [— Q,Q], and K a (t) = (sin Qt)/irt. In order for x(t) to be recovered 
(by an iterative process) they require that / have an inverse over the 
range of x and that is essentially what we require. The inverse of an an- 
alytic / is more complicated, but the fact that the Fourier transforms of 
z(t) and f\z(t)\ vanish over (— <»,0) simplifies the recovery problem. 

Let us write 

w = f(z) (38) 

and 

z = <p(w) (39) 

for the inverse and think first of the problem of recovering z (t ) from w (t ). 
In case / maps \z \ < m one-one onto some region D*, there is no problem 
since <p is single valued over D*. In general tp is not single valued and we 
have to know something about z (r) in order to decide what element of 
<p is "the" inverse. For example, suppose f(z) = 2z + z 2 . Then given 

w{t) = 2ae u + a 2 e i2t 

we do not know whether z(t) = ae lt or z(t) = —2 — ae lt without some 
additional knowledge, such as for example, lim u — „ z{t + iu) = 0, or \z(t) \ 
< 1. In this example the inverse function, 

2 = tp(U)) = -1 + (1 + U)) 1 ' 2 

is not a single-valued function of the complex variable w and one gen- 
erally speaks of two branches of the inverse function. The branches have 
singularities at w = — 1, the image of 2 = — 1 where f'(z) = 0. Clearly in 
this case, if we require \z(t)\ < 1 then we know 



z(t) = -l + Vl+w(t) 
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where VT = 1. Then <p\w(t + iu)} is bounded for u > and analytic for 
u >0. If we relax the requirement to \z(t + iu)\ < 1 forw > b then tf\w(t 
+ iu)\ will be analytic for u > b, but not necessarily for u > 0. 

In general, we require that z(t) and f(z) are so constrained that 
(p\w(t)\ is analytic for u > uq. 

4. 1 Received signal in the range of the inverse function 

Now we are not given w(t) but instead we have a filtered version w a (t). 
Suppose <p\w(t)} is analytic in the upper half-plane u > and w a (t) is 
sufficiently close to w(t) that <p\w a (T)\ is also analytic in the uhp u > 0. 
We say then that w a (t) is in the range of the inverse function. A simple 
sufficient condition for this is that <p be an entire function. Also the 
channel could have sufficiently large bandwidth for w a (t) to be close 
enough to w(t). 

We assume then that 

w a (t + iu)£D*. u>0 (40) 
where 

tp(w) is analytic for w £ D* (41) 

\tff(w)\ ^Mforw £D* (42) 
Then we may take the inverse of w a {t) to obtain 

z a (t) = <p\w a {t)\. (43) 
Now we will see that the Fourier transforms of z a (t) and z{t) agree over 

First, it follows from (35)-(37) and Ex. 1 of Sec. 2.2 that the Fourier 
transforms of w(t) and w a {t) agree over (—<*>,«); i.e., the Fourier trans- 
form of \w(t) — w a (t)\ vanishes over (— °°,a). Then from Theorem 1 we 
have 

\w(t + iu) - w„(t + iu)\ < e~ au sup \w(t) - w a (t)\. (44) 

t 



From (42) we have 

| <p{w) — if(w (t )\ < M\w — w a \ for w £ D* 

w„^D* 



(45) 



Thus [z(t) — z a (t)\ is the boundary value of a function bounded and 
analytic in the uhp satisfying 

\z(t + iu) - z a (t + iu)\ < Me~ au sup \w{t) - w a (t)\. (46) 

t 
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Hence from Theorem 1, 

z(t) - z a (t) = h a (t) (47) 

where the Fourier transform of h a (t) vanishes over (— °°,a). Since the 
Fourier transform of z(t) vanishes outside [0,ft] and a > Q., (47) implies 
that we can bandlimit z a (t) with an appropriate low-pass filter to obtain 
z(t). Thus if Ka ttx (t) is any kernel of Li satisfying 



I 



K Qia {t)e- iwt dt = 1, 0<co<il 

0, co>a (48) 



we have from Ex. 6, Sec. 2.2, with the convention that z{t) is continu- 
ous, 

z{t) = C " z(s)K nia (t - s)ds (49) 

and since the Fourier transform of h a vanishes over (— °°,a) we have (Kn, a 
®h a )(t) = 0;i.e., 

z(t) = C ~ z a (s)K Q<a (t - s)ds. (50) 



4.2 Pre-detection filtering 

In case the received signal w a (t) is not in the range of the inverse 
function, we may under suitable conditions recover z(t) by appropriate 
filtering before (and after) detection. Here then we replace (40) with the 
condition 

w(t + iu) G D* for u > u (> 0). (51) 

It follows from (44) and (42) that for sufficiently large 6 we have 

w a (t + iu) GD[ for u > b (52) 

where D\ is slightly larger than D* and 

(p(w) is analytic for w (E D[ (53) 

\<S(w)\ < Mi for wE.D\. (54) 

Then we have w a (t + ib) in the range of the inverse function. 
Now the Poisson kernel with parameter u 

P u (t)=-—^, u>0 (55) 

reproduces functions bounded and analytic in the uhp from their 
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boundary values. (See the proof of Theorem 1 in Appendix A.) We 
have 

w a (t + ib) = C w a (s)P b (t - s)ds. (56) 

That is, we may determine w„(t) along a line u = b parallel to the real 
axis by convolving w a (t) with the Poisson kernel (parameter u = b). This 
operation we term Poisson filtering. Since, by assumption, w a (t + ib) 
is in the range of the inverse we may take the inverse of w a (t + ib) to 
obtain 

z a (t + ib) = dw a (t + ib)\ (57) 

which is analytic in the uhp and then as argued before 

\z(t + iu + ib) - z a (t + iu + ib)\ < M 2 e~ au . (58) 

So the Fourier transforms of z(t + ib) and z a (t + ib) agree over 

Thus if the conditions (51), (41), and (42) are met we may by suitable 
pre-detection filtering (Poisson filtering) obtain a function z a (t + ib) 
which corresponds to replacing z (t) at the transmitter by z(t + ib); i.e., 
from the reproducing property of the Poisson kernel 

w(t + iu) = f w(s)P u (t - s)ds 

= f\z{t + iu)) = f ~ f\z(s)\P u (t - s)ds. (59) 

Then interchanging the order of convolutions with Pb(t) and K„(t) in 
(35) we have 

w a (t + ib) = C f\z(s + ib)\K a (t - s)ds. (60) 

This relation has been noted by Foschini. 8 

The Poisson kernel is a contraction operator; i.e., it averages the values 
of a function so that the range of the resultant is no larger than the range 
of the function. Also as a filter it has the frequency response 

XPu(t)e- iu "dt = e- u \"\, (u >0). (61) 

We have 

z(t + ib) = f~ z(s)P b (t - s)ds (62) 

and therefore for large b we would expect the range of z(t + ib) to be 
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appreciably less than the range of z(t) for a wide class of z(t) since for 
very large b it is only the low-frequency content of z(t) that contributes 
appreciably to z(t + ib). It is possible, though, for the low-frequency 
content to be such that, for every u > 

z(t + iii) ~ cos y/t + i sin V7 as t — *■ °°. 

If we require z(t + ill) to tend uniformly to a limit z as u — °°; i.e., 

\z(t + iii) - z \ <e(u), -°° <£<«>, (63) 

where c(u) -» as u -* °° and in addition 

/'(2o) * 0, (64) 

then p(ic) will be analytic in the neighborhood of w = f(z ) and so (51) 
will be satisfied for sufficiently large u with (41) and (42) holding. 

The condition (63) is not a severe constraint. In fact, all the simple 
sufficient conditions, given in (21)-(29), for x(t) to have a Hilbert 
transform imply (63) with z = 0. However, x(t) may have a Hilbert 
transform without (63) holding. 

In connection with pre-detection filtering, we note that equivalent 
Poisson filtering can be effected at the carrier frequency (or an inter- 
mediate frequency) in the receiver before the synchronous demodulation 
of the received signal o R (t) indicated in (34). That is, if the signal o R (t) 
is passed through a filter whose frequency response Fb (to) satisfies 

F b {oi) = e-*<"-"»>, co > co c (65) 

F 6 (-co) = F 6 (co) 
a signal OR(t;b) is obtained which we may identify as 

a R (t;b) = Re [e^wjt + ib)\. (66) 

Then synchronous in-phase and quadrature detection of <tr(£;6) yields 
the real and imaginary parts of w a (t + ib). Thus the Poisson filtering 
may be accomplished with a single equivalent frequency-translated filter, 
whereas the direct Poisson filtering of the complex signal w a (t) requires 
two Poisson filters acting separately on the real and imaginary parts. 

4.3 Post-detection filtering 

When Poisson filtering is required to bring the received analytic signal 
within the range of the inverse function the low pass filtering after de- 
tection must be modified. The output of the detector is z a (t + ib) and 
we need z(t). Now the Fourier transforms of z(t + ib) and z a (t + ib) 
agree over (— °°,«) and 

z(t + ib) = f " z(s)P b (t - s)ds. (67) 
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The Poisson filtering operation has an inverse for bandlimited func- 
tions; i.e., 

z(t)= f° z(s + ib)Q b (t - s)ds (68) 

where Qb(t) is any function of Li satisfying 

f Q b (Oe-' a,£ di = e fe <", < co < £2. (69) 

Since the Fourier transforms of z(t + ib) and z a (t + ib) agree over 
(— °°,a), the Fourier transform of 

z(t) ~ f z a (s + ib)Q b (t - s)ds (70) 

vanishes over (—<»,«), and since the Fourier transform of z(t) vanishes 
outside [0,0], we have [cf. (47)-(50)] 

z(t)= C~ z a {s + ib)k(t-s)ds (71) 

where k(t) = k{t;b,$l,a) is any kernel in Li satisfying 

k(t)e- iut dt = e" b , < a) < fi (72) 



s. 



= 0, co > a. 



That is, the post-detection filtering must invert the pre-detection fil- 
tering over the band [0,fl] and remove frequencies greater than a. (Of 
course, in a practical system we are interested in recovering only x (t) so 
that only one post-detection filter is required, acting on the real part of 
z a (t + ib).) 

V. SPECIALIZATION TO BAND-PASS SIGNALS 

In case the base-band signal x(t) is of the band-pass type, i.e., a signal 
whose Fourier transform vanishes outside [rQ,Q] and [— !),— rft] where 
< r < 1, the detection theory may be modified so that no Poisson fil- 
tering is required. In this case the inverse function may be replaced by 
an entire function, in particular, a polynomial. All we require for the 
recovery of band-pass signals is that /|z(t)} be analytic in the uhp and 

/'(O) ^ 0. (73) 

Then 

w = f(z) = Y, a k zk f° r | z | sufficiently small (74) 

k=o 
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2 = <p(w) = E bk(w — ao) k for \w — ao\ sufficiently small. (75) 

Hflnfci-(ai)- 1 -|f(0)]- 1 . 

Now for band-pass signals x(t), the Fourier transform of the analytic 
signal 2 (t) vanishes over (— «>,rft) where Q > and < r < 1. Thus by the 
Paley-Wiener Theorem for L m , 

\z(t + iu)\ <e- r «"sup \z(t)\, u>0. (76) 

t 

Then we have 

w(t + iu) — ao = Y, a k\z{t + iu)\ k for sufficiently large u. (77) 

k=\ 

It follows that the Fourier transform of \w(t) — a ] also vanishes over 
(— oo,rfi). Hence the Fourier transform of \w a (t) — o,q\ vanishes over 
(-°°,rfl). 

Now let <p*(w) be any entire function of the form 

<p*{w)= f c k (w-a ) k (78) 

/e=l 

where 

c k = b k , for /z = 1,2, . . . ,n (79) 

and rc is an integer such that 

nr > 1. (80) 

Defining 

z* a (t) = <p*\w a (t)\ (81) 

we have 2*(t) analytic in the uhp and for sufficiently large u 

z(t + iu) — z* a (t + iu) 

= (p\w(t + iu)} — (p*\w a (t + iu)\ 

= E b k [\w(t + iu) - a \ k - \w a (t + iu) - a } k ] 
k=i 

+ E b k \w{t + iu) - a \ k 

k=n+l 

- E c k {w a (t + iu) - a \ k . (82) 

k=n+l 

Since the Fourier transforms of \w(t) - a \ and \w a (t) - a \ vanish over 
(-oo,rfl), the last two sums in (82) are of the order of exp \-(n + l)rfiu}. 
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Also, for k > 1 

\\w{t + iu) — ao\ k — \w a (t + iu) — ao\ k \ 

= 0\\w(t + iu) - w a (t + iu)\\ (83) 

= 0(e~ au ). 

So it follows from the Paley- Wiener Theorem (with Theorem 3) that 
the Fourier transforms of z(t) and z* a (t) agree over (— °°,B) where 

B = min \a,(n + l)rfi} > fi. (84) 

Therefore 

z(t) = f~ z *MK QiB (t - s)ds (85) 

where Kqb is any function of L\ satisfying 

Koji(t)e- iut dt = 1, 0<a)<fi (86) 



x: 



= 0, a> > B. 



Thus for band-pass signals we have the option of replacing the inverse 
function <p(w) by an equivalent entire function <p*(w) so that it does not 
matter whether or not the received analytic signal w a (t) is in the range 
of the inverse function <p(w). In particular, <p*(w) may be a polynomial 
of degree n where n is roughly the ratio of the upper and lower cut-off 
frequencies of the base-band signal. 

VI. DETECTION OF EXPONENTIAL MODULATION 

The exponential modulation law/U) = e z offers the unique advantage 
of eliminating the need for preliminary in-phase and quadrature de- 
tection of the received single-sideband signal OR(t). In this case we 
have 

z = <p(w) = log w (87) 

or using Log to denote the real part of the logarithm, 

x(t) - Log \w(t)\ (88) 

y(0 = argMi)I. (89) 

We may regard either x(t) or y(t) as the signal to be recovered. 
The transmitted signal is 

a(t) = Reexp\iw c t + z(t)\ (90) 

and 

a(t) + ia(t) - exp [ia c t + z(t)} (91) 
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where a(t) denotes the Hilbert transform of o{t). The envelope of a(t) 
is 

E\a(t)\= \a(t) + ia(t)\ = e* (t) (92) 

The instantaneous phase of a(t) is 

A\a(t)\ = arg \<x(t) + ia(t)\ = u> c t + y(t) (93) 

For perfect transmission; i.e., <r/?(t) = o(t), we have 

x(t) = Log E[a(t)} (94) 

and using an ideal discriminator (FM detector) we obtain 

y'(t)=^-A\o(t))-u c . (95) 

at 

Replacing a(t) by orH) we have 

x a (t) = Log E\a R (t)) = Log \w a (t)\ (96) 

y^t) = -f A\a R (t)\ - a> c = — Arg |w a (t)}. (97) 

at at 

Now if w a (T) is zero-free in the uhp, then 

z a (r) = log w„(r) (98) 

is analytic in the uhp and by the previous theory the Fourier transforms 
of z(t) and z a (t) agree over (-«>,«). In this case the Fourier transforms 
of x a (t) and x(t) agree over (-a,a). Also the Fourier transforms of y' a {t) 
and y'(t) agree over (-a,a). So if the received analytic signal w a (t) is 
zero-free in the uhp, x(t) may be recovered by taking the Log of the en- 
velope of the received signal and then filtering with an ideal low-pass 
filter having unity transmission in the band [— fl,Q] and zero transmission 
outside the band [-a,a\. Similarly y'(t) may be recovered by filtering 
the output of an ideal discriminator acting on the received signal. 

Later, in examining the bandwidth requirements of single-sideband 
exponential modulation we give sufficient conditions for w„{t) to be 
zero-free in the uhp so that the simple detectors described above may 
be used. 

The simple detectors can always be used with appropriate pre-de- 

tection and post-detection filters, since for sufficiently large b, w a (,r + 

ib) will be zero-free in the uhp. We can give an estimate for 6 under the 

condition 

sup|2(t)|<m. (99) 

t 

We have 

e~ m < \w(t)\ <e m (100) 
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and consequently 

e~ m < \w(t + iu)\ < e m , u > 0. (101) 

Also, w(t) is zero-free in the uhp and hence [cf. (219)] if 

\w(t + ib) - w a (t + ib)\ < e~ m (102) 

then w„(t + ib) will be zero-free in the uhp. From (44) we have 

\w(t + ib) - w a {t + ib)\ < e~ ab sup \w(t) - w a (t)\ (103) 

t 

and since 

w a (t)= f w(s)K a (t - s)ds (104) 

we have 

sup \w a (t)\ < sup \w(t)\ • \\K„\\i < e^WKJi (105) 

t t 

and hence 

sup \w(t) - w a (t)\ < {1 + ||K„|| 1 }e m . (106) 

Thus if 

e-« b \l+ ||/f a || 1 }<c-*» (107) 

then (102) will be satisfied. That is, w„{t + ib) will be zero-free in the 
uhp for 

, v 2m + log|l + \\KJA 

b > ^ } -^ m (108) 

a 

Thus if (99) is satisfied we may use a frequency-translated Poisson filter 
with parameter b satisfying (108) to obtain a R (t;b) [cf. (66)] and then 
operate on the envelope and phase of a R (t\b) as before. Then the ap- 
propriate post-detection filtering may be employed to recover x(t) and 

y'(t). 

VII. NOTE ON THE FACTORIZATION OF CERTAIN POSITIVE FUNCTIONS 

The detection theory for SSBEM has important application to the 
problem of factoring certain positive functions of exponential type, i.e., 
certain positive bandlimited functions. 

Voelcker 22 has proposed a scheme for demodulating conventional 
single-sideband signals via envelope detection. Conventional SSBAM is 
characterized by linear modulation; i.e., f(z) = z. There is no bandwidth 
expansion so we assume that the Fourier transform of z(t) vanishes 
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outside [0,fi] and the channel is such that the received signal is simply 
the transmitted signal; i.e., 

a R (t) = a(t) = Re e^zit). (109) 

The envelope of the received signal is \z{t)\ . Voelcker's scheme requires 
first that z (r) be zero-free in the uhp in order that z (t ) may be recovered 
from \z(t)\. This is insured by requiring Re z(t) = x(t)> 0. We have 

z(t + iu) = f " z(s)P u (t - s)ds (110) 

and since the Poisson kernel is positive, 

Rez(t + iu)= f" x(s)P u (t - s)ds > 0. (Ill) 

Therefore 2 (t) is zero-free in the uhp. The function log z(t) is analytic 
in the uhp and with some additional conditions on z(t), e.g., 

\imz(t + iu) = 1, (112) 

the imaginary part of log z(t) can be determined from Log \z(t)\ and 
hence z(t) can be recovered from \z(t)\. In particular, if 

x(t) = l+g(t), (113) 

where 

g(t) > -1 and g(t) belongs to L p (1 < p < °°), (114) 

then Log \z(t)\ will belong to L p and will therefore have a Hilbert 
transform. A more attractive condition for recovering x(t) is the condi- 
tion 

g(t) > -1 and g(t) of band-pass type. (115) 

Then if (115) is satisfied, the Fourier transform of \z(t) - 1} vanishes 
outside [ri2,H], where < r < 1, and hence 

w{t) = log z(t) (116) 

is analytic in the uhp and satisfies 

w(t + iu) = log [1 + \z(t + iu) - 1}] (11?) 

= 0\\z(t + iu) - 1|| = 0(e- rttu ), u — °o. 

Therefore, if (115) is satisfied, the Fourier transform of w(t) vanishes 
over (— °°,ri2) and hence the Fourier transforms of Log \z(t)\ and 
arg \z(t)\ vanish over (-rfl.rfl). That is, if (115) is satisfied, then the log 
of the envelope of z(t) and the phase (arg) ofz(t) are high-pass func- 
tions. 
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If either (114) or (115) are satisifed we have 

arg| 2 < t )l^(0=-r^^U (118) 

and then 

2(0= |«(0|eH**>. (119) 

The practical problem encountered here is in approximating the 
Hilbert transform in (118). The function Log \z(t)\ is not band-limited* 
and the implementation of (118) requires a filter whose frequency 
characteristic is, formally, 

H(u) = I — — dt = -i I dt = -i sgn co. (120) 

•/-<» irt J-*, -nt 

( — 00 < w < 00 ) 

These stringent filter requirements can be avoided, for we can, by proper 
application of the previous theory, ignore the frequency content of 
Log | z (0| outside the band (-«,«) where a > fi. Actually if (114) is 
satisfied we may take a = ft. 

The Hilbert transform problem is simplified if we begin with a filtered 
version of Log \z(t)\; viz., 

MO = J Log \z(s)\h a (t - s)ds (121) 

where h„(t) is an even real-valued function whose Fourier transform 
satisfies 



fia(w) = f h a (t)e- iwt dt =1 for -a < a> < a. 



(122) 



We suppose further that /i„(0 is sufficiently smooth to have a Hilbert 
transform h a {t). Then the Hilbert transform of \ a (t) is given by 

X«(0 = f> a (t) = C Log \z(s)\h a (t - s)ds (123) 

= I <p(s)h„(t — s)ds. 

Then defining 

w„(t) = \ a (t) + itp a (t) (124) 

we have 

w a (t) = C ^ log \z(s)\h a (t s)ds. (125) 



f Unless z(t) = constant. See Theorem 6 in Section IX for a stronger statement. 
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Then according to the previous theory 

z(t) = f " k(t-s) exp \w a (s)\ds (126) 

where k(t) is any function of Li satisfying 

k(t)e~ iut dt = 1, 0<co<fl (127) 



x: 



= 0, oj > a. 



In case (115) is satisfied the function h a {t) need satisfy (122) only over 
the intervals (rtt,a) and (-a,-rfi), since the Fourier transform of 
Log \z{t)\ vanishes over (-rfi,rfi). That is, h„(t) can then be chosen so 
that the equivalent Hilbert transform kernel h„(t) has a Fourier trans- 
form that is more easily approximated (within a linear phase factor 
e~ io)T ) by practical filters. 

VIII. NOTE ON LOGARITHMIC COMPANDING 

Suppose g(t) is a function belonging to L p for some p satisfying 1 < 
p < <» and suppose the Fourier transform of g(t) vanishes outside [—12,12]. 
Companding functions / are sometimes used to compress the range of 
g{t) for transmission; i.e. f\g(t)\ is transmitted rather than g(t). Landau 
and Miranker 11 - 12 showed that g(t) (in L 2 ) can be recovered from the 
bandlimited version of f\g(t)\ with suitable conditions on /. The recovery 
is accomplished by an iterative scheme. Here we use the detection theory 
to give an explicit solution to the problem of Landau and Miranker for 
the case 

/(*)=- Log (1 + x), x>-l. (128) 

Accordingly, we further require g(t) to satisfy 

g(t) > -1. (129) 

The function f(x) given by (128) is not an odd function, as one might 
desire for companding purposes, but is interesting because the recovery 
problem is simple. 

The fact that g(t) is a bandlimited function belonging to L p for some 
p satisfying 1 < p < <» implies [cf. (29)] 

m- r *) ***-* * (180) 

J-*- ir(t — s) 

from which one can conclude with the aid of Holder's inequality that 

lim g(t) = 0. (131) 

t— ±eo 
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It follows from (129) and (131) and (128) that 

f \g(t)\ belongs to L p . (132) 

Since g(t) is bounded and belongs to L p it follows that g(t) belongs to 
L P ' for every p' satisfying p < p' < °°. Hence f\g(t)\ also belongs to L p < 
for such p'. 

Now we suppose we are given 

M*>- f'/te^l ^f-'^ (133) 

*/-°° 7r(t — s) 

where the integral is absolutely convergent by Holder's inequality. In 
fact (Ref. 16) Xa(t) belongs to L p , and therefore has a Hilbert transform. 
Furthermore, since \n(t) is bandlimited, its Hilbert transform is given 
by 

X fl (t) = *a(*)= rxds ^'^^'^ ds (134) 

J— Tr(t — s) 



I 



i - cos n(t - s) 

f\g(s)\ ds. 

ir(t — s) 



Defining 



we have 



where 



w n (t) = \ a (t) + i<p n (t) (135) 



w Q (t) = C " f[g(s)\K u (t - s)ds (136) 



pint _ i 
K a (t)= , . (137) 

I7r£ 



So m)j2(t) is an entire function which is bounded in the uhp. 

Now (1 + g(t)\ is a positive bandlimited function which can be repre- 
sented as (Theorem 7.5.1 with Theorem 6.4.5, Ref. 6) 

l+g(t) = y(t)y(t) (138) 

where the Fourier transform of y(t) vanishes outside [— fi/2,fi/2] and 7(7-) 
is zero-free in the uhp. Then z(t) defined by 

z(t) = y(t)e iilt/2 (139) 

is a function whose Fourier transform vanishes outside [0,ft] and z(t) 
is zero-free in the uhp. Thus we have 

1 + g(t) = \z(t)\*. (140) 
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We may assume that [cf. (131)] 

lim z(t) = 1. (141) 

Then z(t) is given by 

z(t) = exp w(t) (142) 



where 



and 



ds 



'(r)=^ f°Log{l+g(s)} - ' - (143) 

2 •/-« 7r(r — s) 



lim w(£ + iu) = \(t) + i>(t) (144) 

u—G+ 

X(*) = | Log {1 +*(*)) (145) 

^(t) = \(t), the Hilbert transform of \(t). (146) 

We see from (136) and (143) that 

w Q (t) = C " u;(s)Ko(t - s)ds (147) 

= J^log{z( S )}K n (t-s)ds. 

Then the Fourier transform of zn(fc) defined by 

z (t) - exp {i0 O (t)) (148) 

agrees over (-°°,Q) with the Fourier transform of z(t). Since \z(t) — 1} 
belongs to L p and its Fourier transform vanishes outside [0,ft] we 
have 

.(0-1- f"l«oW-l|^^«fc (149) 

J-co 7r(t — s) 

Writing z (t) = x(t) + iy(fc) we have 

x(t) = 1 + f " {e x °<*> cos «,(s) - 1} sm ° ( *~ g) <fe (150) 
J-a= 7r(t — s) 

yit )= f%x«( S ) S m^( g ) sin " a " s) d g (151) 

Then we have 

*(t) = * 2 (t) + y 2 (*)-i. (152) 
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% 



Thus the recovery problem is solved by means of the Hilbert transform 
in (134) and the formulas (150)-(152). 

Expressing the solution in terms of the bandlimiting operator j8n, 
defined for h in L p , 1 < p < °° , by* 

. , x C , , s sin Sl(t — s) , 

Mt)= h(s) -f -^ds, (153) 

*/-» 7r(t — s) 

the resulting class of functions denoted by B p (fi), we have 

g(t) = \% n exp {X fi (0 + i"X fi (0}| 2 - 1 (154) 

where X^it) is the given function 

Xn{t)=ls Q Log\l+g(t)\ (155) 

ginB p (fi), g>-l, 

and Xn is the Hilbert transform of Xq. 

The solution (154) is deceptive in that it suggests that Xn(t) may be 
any function of B p (fi), 1 < p < co , since g(t) given by (154) is a function 
of Bp(fl) satisfying g > —1. However the solution was obtained on the 
premise that X n (£) is a given function of the form (155). All functions in 
B P (Q) do not have the representation (155). The crucial point is that the 
function 

z(t) = Snexp\Xdt) + i\n(t)\, \z(t) - l}inB p (fl), (156) 

whose Fourier transform vanishes outside [0,fl] should extend as a 
function zero-free in the upper half-plane. Then, and only then, ac- 
cording to the general theory, will we have 

J8 n log \z(t)\ = X Q (t) + iX a (t) (157) 

and hence 

% a Log\z(t)\=^Z a Log\l + g(t)\ = X u (t), g in B p (fi), (158) 

On the other hand, if (158) is known to hold, implying (157), then z (t) 
must necessarily extend as a function zero-free in the upper half- 
plane. 

We state this important result as 

Theorem 4. Given a function X n (t) in B P (Q), for some p satisfying 1 < 



f The operator % a can be extended to certain other classes of functions. For example, Sa 
is an identity for the constant function, which fact is used in (154). 
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p < », the equation (155) has a solution g(t) in the same class B P (Q) 
satisfying 

g(t) > -1 
if and only if the function 

*(t)-£o«cp{Xo(t) + rto(t)) 

where \n(t) is the Hilbert transform of \q(£), extends as a function 
zero-free in the upper half-plane. Then the solution of (155) is given 
by (154). 

IX. BANDWIDTH REQUIREMENTS FOR EXPONENTIAL MODULATION 

We have seen that for a wide class of analytic signals z (t) and modu- 
lation laws f(z) the bandwidth requirement for transmitting f\z(t)\ and 
recovering z (t) is ft + c (for any e > 0) where ft is the bandwidth of z(t), 
provided we allow the use of Poisson filtering at the receiver. In case the 
inverse function z = <p(w) is an entire function there is no need for 
Poisson filtering. If we look at the overall system design, as contrasted 
to a detection problem, it is reasonable to ask for the bandwidth re- 
quirements for a given f(z) and a fixed receiver, namely the inverse 
function <p(w) followed by a low-pass filter, such that we recover all z(t) 
whose Fourier transforms vanish outside [0,ft] and satisfy some sort of 
norm constraint, say \z(t)\ < m. The problem then is to specify a channel 
of finite bandwidth, i.e., a function K a< p(t) in L\, with ft < a < /3, satis- 
fying 

R a ,p(^) = j^K a ^(t)e-^dt=l, 0<o><« (159) 

= 0, w > 
such that the received (bandlimited) analytic signal w a< p(t), given by 

w a , p (t) = f " f\z(s))K a>p (t - s)ds, (160) 

satisfies 

<Awa,p(T)} analytic in the uhp (161) 

for all z(t) whose Fourier transforms vanish outside [0,ft] and which 
satisfy 

\z(t)\<m, -co<t<<*>. (162) 

We assume of course that f(z) is analytic for \z\ < m. We would like to 
make the channel bandwidth /3 as small as possible consistent with (161) 
and (162). Clearly, we may take ft = 1 with no loss in generality. We de- 
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fine the minimum bandwidth Po(m) as 

0o(m) = inf0(m) (163) 

where the infimum is over all functions K a #(t) subject to (159), (161), 
and (162) with fi = 1. In taking the infimum we may allow a — 1. 

The determination of /3 (m) is in general a very difficult problem. We 
give some estimates here for /3 (m) for the case/(z) = e z , <p(w) = log w. 
In this case, (161) is satisfied if and only if 

Wa,/3(-r) is zero-free in the uhp. (164) 

We have 

w(t) = exp\z(t)} , (165) 

and with (162) 

e~ m < \w(t)\ <e m . (166) 

If w ai ff(t) is sufficiently close to w(t), (164) will be satisfied; i.e., a 
sufficient condition for (164) is 

\w{t) - w a , p {t)\ <e~ m . (167) 

It is intuitively obvious, with the freedom we have in defining R. a ,p, 
that for sufficiently large /3 we can find a function K a ^(t) such that 
w ai p(t) given by (160), with /(z) = e z , will satisfy (167). It is important 
to note in this connection that, although the definition of i£ a>/3 (co) for a> 
< does not affect w ai p{t), we are free to define it ai/3 (a>) for a> < (as well 
as for a < to < /?) in the most favorable way to obtain the estimate 
(167). 

First we obtain lower bounds for Po(m)- 

9.1 Lower bounds for 3 (m) 

We can obtain a lower bound for /?o(w) by taking 

z(t) = me". (168) 

We have 

00 m^ 
w(t) = exp {me 11 ) = £ —e ikt . (169) 

k=o k\ 

Now assume the channel cutoff frequency /3 satisfies 

n < (3 < n + 1 (170) 

where n > 1 is an integer. Then 

w a At) = 1 + me lt + E a k e ikt (171) 

k = 2 

ANALYTIC MODULATION SYSTEMS 523 



where the a* depend on m and the definition of R a ^(o)) for 1 < a> < n. 
We have 

w a ,0{t)= ft (l-f) (172) 

k=i \ h' 

where f = e lt and 

E—=-m. (173) 

We require w ai p(t) to be zero-free in the uhp; i.e., 

|f*|>l. (174) 

Thus 



™ = £ — 



* t 777 ^ »• (175) 



Therefore, we must have /3 > n > m in order for w a ^(t) to be zero-free 
in the uhp. Then 

/3 (m) > M+ (176) 

where [m] + is the smallest integer which is not less than m. 

9.2 Lower bound for small m 

We know that Po(m) > 1 for any m > but (176) does not say how 
much @o(m) must exceed 1 for < m < 1. For sufficiently small e and 
correspondingly small m we can show that 0o(m) > 1 + e. 

For small m, we have 

z 2 (t) 
w(t) = 1 + z(t) + -^ + 0(m 3 ). (177) 

Now by Corollary 2 of Theorem 2 the Fourier transform of z 2 {t) vanishes 
outside [0,2] and we would like to find a z(t) such that a channel filter 
with a sharp cut-off, i.e., (3 = 1 + e, acting on a small z 2 (£) gives a large 
(negative) output at t = 0. For sufficiently small c and fixed m we can 
accomplish this by taking 

m pit/2 

z(t) = ^{l + iSAet)} (178) 

v2 

where S n (t) is a sine polynomial, 

Sn(t) = £ a* sin fc£ (179) 

A=l 
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with real coefficients a* and 

max |S n (t)| =1. (180) 

Also we require 

ne = - (181) 

2 

so that the Fourier transform of z(t) given by (178) vanishes outside [0,1]. 
Also we have 

max \z(t)\ =m (182) 



and 



m pit 12 

w{t) = l--^-\l + iS n (et)) 



m< _. 



where 



+ — e*|l -Sfret) + 2iS n (et)} (183) 

4 

+ Rs(t), 



- z k (t) 
Rait) = L — y (184) 

k = 3 kl 

m 3 1 
\R 3 (t)\< — , (m<4). (185) 



(-7) 



We have 



J" me"' z 

_ wtt)if^(* - $)<** = 1 — ^=- {1 + S„(«t)} 

^j2 y-» CO 

+ — J_ m e'*{l ~ S'fcfl + 2iS n (e£)}K af/3 (t - Od£ (186) 

and w a> ff(t) is a polynomial of degree In in exp (iet) where e = l/2rc. We 
have 

2ie lt S n (et) = -e lt £ o*e-'** d (187) 

+ e u f; a A e* feet 

2n 

e ft Sj(«0 = e" L foe 1 '*". (188) 

/e=-2n 
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So 



w a ,p(t) = 1 - 



me 



it/2 



V2 



\l + iS n (et)} 



4 



where 



1 - e lt £ b k e iket - e lt £ a k e- iktt 

k=-2n k = l 

+ r 3 (t) 



r 3 (t) = f^Rd&Kajit ~ Sm 



(189) 



(190) 



Since Rz(t) is a periodic function of the form Sq c^e ik(t we may take 
for K ai p(t) any function of Li whose Fourier transform satisfies 

RaA") = 1> < CO < 1 

= 0, co > = 1 + e ' 



(191) 



It is shown in Appendix B that there exists a function K a< p(t) whose 
Fourier transform satisfies (191) with 

§~ \K a , p {t)\dt < 1 + i log (l 4- ^) • (192) 



Thus 



|r 3 (t)| < max \R 3 (t)\ C " |X„,^(i)|dt 



! (-7) 



We have from (189) 



... . m m' 
w atfi {Q>) = !- — + — 



1 - L bft " L a* 

fe = -2n fe=l 



+ r 3 (0). 
Now S„(£) is an even function and S n (0) = 0. Thus 

2n 



(193) 



(194) 



E &fc = 0, b-k = bk 

k=-2n 


(195) 


-1 2n _fc n 
k=-2n k=l 2, 


(196) 


£ h - b ° 
k=-2n 2 


(197) 
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and 



b = ^- Cs*(t)dt = l£al<l. (198) 

2ir J-* 2 k=i 

Now we may take S n [t) to be an approximation to sgnjsin t}. In particular 
we may take S n (t) to be the function found by Szego 19 which maximizes 
2" ak subject to (180). He gives the inequality 

Za k <M n = -^- £ cot( kl [ )~-logn, n - «. (199) 
i n + 1 i<fc<n \2(n + 1)/ 7r 



A odd 



Here we should not identify ak with the terms in the second sum. For 
equality in (199) we must have 



S n (-^-) - 1. l<k<n, k odd. 
\n + 1/ 

Equality in (199) is attained for 



(200) 



(201) 



k odd 



where 



. n + 1 



/„(*) = 



sin 



(202) 



(n + 1) sin - 



It is easy to show that 

£ I n It- -^-) =1 for n odd 
J= _ n \ n + 1/ 



k-- 

(k odd) 



* = -(n-l) \ fl + 1/ 2 



(203) 



(204) 



(A odd) 



Hi 



for n even 



It follows from (203), (204), and (201) that S n (t) given by (201) satis- 
fies 

-l<S n (t)<l. (205) 
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From (201) and (202) we find that 

sin z - 
4 / k \ 2 nodd 



. kir 
sin z — 



a k 



•—(i-—) 

n + l\ n + l/ 



. kit k = l,2,...,n (206) 

sin 

n + l 

(nkir \ 2 
sin— ) 

\ n GVen (207^ 

°* n + l\ n + l) . fe x ' & - 1,2, . . . ,n 

sin 

n + l 

It is shown in Appendix C that M n given by (199) satisfies 

M n > - log n + - Aog - + 7) (208) 

7T 7T \ 7T / 

where 

T = 0.5772. . . (Euler's constant). (209) 

So with S n (t) given by (201) and z(t) given by (178) we have from 
(193), (194), (197), and (208) 

/n\ ^ 1 m , m2 [-, b 2, 2/, 4 \ 

^(0)<l-- + T { 1 - T --lo gn --(log ; + 7 )} 

+ — =*— {l + ik,(l + *)] (210) 

3!(x-2)' * V 3/( 

where 1 < a < 0, 

c 2n 



1 C* 

bo = — I S*(£)ctt ~ 1 as n 

27T \) —K 



Now if we set 

it is clear that 
Since 



in* 

— logn = l (211) 

Z7T 



w aij3 (0) < for sufficiently large n. (212) 



2rc 

Wa,/?(t) = Z c k e ikt ' 2n (213) 

fe=0 
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where the Ck are real and c > 0, we have w a ^(iu) real and 

lim w ai0 (iu) = Co > 0. (214) 

So (212) and (214) imply that w at0 (t) has at least one zero on the positive 
imaginary axis; i.e. for sufficiently large n, and m given by (211), we must 
have > 1 + l/2n (to pick up at least another harmonic) in order for 
w a ,p(t) to be zero free in the uhp. Thus 

0o ( V -^-) > 1+ — for sufficiently large n. (215) 

\ log nl 2n 

In connection with obtaining lower bounds for /3 the idea comes to 
mind that we might be able to find a z (t ) satisfying |z(£)| < m (for suf- 
ficiently large m) such that the Fourier transform of w(t) would vanish 
over a large interval (1,0). Then if w a ,p(t) were not zero free in the uhp 
we would have 0o > 0(m). The idea is to obtain a w ai p(t) which would 
be independent of the choice of K„ t p(t). However, we cannot make the 
Fourier transform of w(t) vanish over large intervals unless z(t) = con- 
stant. 

Theorem 5. Suppose z(t) is a bounded continuous function whose 
Fourier transform vanishes outside [0,0], and suppose that the Fourier 
transform ofw(t), where 

w(t) = exp \z(t)\, 

vanishes over (a,b) where 

a > and b — a > ft. 

Then 

z(t) = constant. 

A similar result holds for the logarithmic function. 

Theorem 6. Suppose z(t) is a bounded continuous function whose 
Fourier transform vanishes outside [0,J2] and suppose its analytic 
continuation z(t) satisfies 

\z(t + iu)\ >€>0 foru>0 

-co < t < co . 

Suppose further that the Fourier transform ofw(t), where 

w{t) = log \z(t)\ 
vanishes over (a,b) where 

a > 0, b - a > Q. 
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Then 

z(t) = constant. 
As an application of the last theorem, we may take ft = n and 

Z (t) = ft (i - Xfcc 1 '**) 
fe=i 

Then, assuming | A& | < 1, we have 

log z(t) = - L-^— 

m=l tfl 

where 

Mm = t (X/e) m 

Then the following is true. 

Corollary. If [\k], k = 1,2, . . . ,n, is any set of n complex numbers and 

f; (X k ) m = for m = p, p + l,p + 2, . . . p + n - 1 

fe=i 

where p is a positive integer, then \k = 0, k = 1,2, . . . ,n. 

Proofs of Theorems 5 and 6 are given in Appendices D and E. 

9.3 Upper bound for fi (m) 

We have 

w(t) = e*M and \z(t)\ < m (216) 

so [wit)] -1 is bounded and analytic in the uhp. Thus the quotient 

wit) w(t) {21V 

is bounded and analytic in the uhp, and is reproduced by the Poisson 
kernel from its values on the real line. Then if 



wit) -w aj p(t) 



wit) 



< 1 for -oo < t < co (218) 



the function w a> pit) is necessarily zero-free in the uhp. Thus a sufficient 
condition for w a ^it) to be zero-free in the uhp is 

\w(t) - w a , p (t)\ <e~ m . (219) 

For lack of something better we will use this condition to obtain an upper 
bound for fioim). To meet this condition for large m will require con- 
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siderably more bandwidth than the lower bound for Po(m). In fact for 
the case z(t) = me lt we have 

wW-Ett**' < 220 > 

k = kl 
and 

w a ,p(t) = 1 + me u + £ a ft e f *' (221) 

2 

where the a* depend on the definition of R a fi{u>) for co > 1. Since 

— = — C \w(t)-w af) (t)\e- ikt dt, k>(3 (222) 

k\ 2tt J-TT 

we have 

YT ^ max | a;(0 - w a ,fl(0| where k > 0. (223) 

Hence in order to satisfy (219) for the case z(t) = me 1 * we must have 

mk ]<e- m . (224) 



max 

fe>/3 



For large m we must have /3 large since 

fe! < V^k k k e- k e l ' l2k . (225) 

We find for large m that 

j8 > pm + o(m) (226) 

where p is the root of 

p = e i+i/p = 3.591121477. 

Thus for large m the upper bound we obtain for Po(m) from the condition 
(219) must be something like 3.6 times as large as the lower bound (~m) 
we obtained previously. We can in fact obtain an upper bound for j3 (m) 
that is close to pm for large m and, as it turns out, is close to the lower 
bound for (io(m) as m -*■ 0. 
To do this we suppose that 

a = n < (227) 

where n is a positive integer. We take 

R a ,e(u>) = 1, < co < n 

= ^^ , n < co < 13 (228) 

|8 — n 

= 0, co > p. 
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^a,/3(w) is defined for a> < in such a way (see Appendix B) that 

An 



r \K a>0 (t)\dt <-\o g 

J-* 7T 



1 + 



Now we write 



where 



3(0 -n) 

w(t) = e^ = P n \z(t)\+ ± ^* 

k=n+l «! 



k=0 k\ 



+ 1. 



(229) 



(230) 



(231) 



Since by Corollary 2 the Fourier transform of P n \z(t)\ vanishes outside 
[0,n], we have 



w 



where 



a,p(t) = J w(s)K a>f }(t - s)d& 
= P n {z(t)) + R n+1 (t) 



R n+ i(t)= C" \ t ^f-\K aif} (t-s)d& 

•/ — [k=n+l Kl J 

\Rn+l(t)\<\ £ Trl f" \K'a.p(t)\dt. 
ik=n+l «! J J-°> 



(232) 



(233) 
(234) 



Thus 

\w(t)-w a ,(,(t)\ < 
Therefore if 



lk=n+l k\ 



1 + 



§_l\K«At)\dt}' (235) 



\K a , (t)\dt < -=Q n (m) 



k=n+l kl 



(236) 



w a ,p(t) will be zero-free in the uhp. Then from (229) and (236) we can 
get an upper bound on /3 for each choice of n. Since 



f " K ai$ (t)e- iut dt =1 for < a < a 



(237) 



it follows that 



J]" \K a , p {t)\dt>l. (238) 
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So the inequality (236) can hold only if m and n are such that the 
right-hand member of (236) is greater than 2. Then setting 

A„(ro) = max (0,[Q n (m) - 2]j (239) 

we see from (229) and (236) that 



-log 

IT 



1+ 4 " 



3K (m ,-,r A " (m) (240) 

defines for each positive integer n an upper bound for Point), viz., 

4 
— n 

<>«<»•> ■ » + e ,JL-i - (24i) 

We take /3„(m) = °° for A n (m) = 0. For fixed m, we have A n (m) > for 
sufficiently large n. Then 

/3 (m) <B(m) = min/3 n (m) n = 1,2,3, .. . (242) 

n 

Since n (m) > n, it is clear that the minimum in (242) will be 0i(m) 
for sufficiently small m. We have 

Ai(m) = 2m- 2 + 0(m- 1 ), m — 0. (243) 

So the upper bound /3i(m) for small m compares favorably with the lower 
bound (215). At least we have the dominant exponential behavior pinned 
down as m — *■ 0; i.e., 

lim m 2 log \p (m) - 1} = -2tt. (244) 

m— 

The function A n (m) defined in (239) behaves like (n + l)!/m n+1 as m 
-*■ and decreases to zero at m = m n where 



(n+1) , 1 , . 1 

m„ m + log (w + 1) + —— 

P 2(1 + p) l + p 



- log 2t + log 



2 2p J 

(245) 

and p = 3.591121477 is defined in (226), 

- = .278464543 
P 

1 



= .1089058525 



1 + P 



2(1 + p) 

-log 2tt + log^^l = -.0219080253. 
2 s 2p J 
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Fig. 1 — Upper bound for transmission bandwidth required in ESSB for simple detection 

of signals z (t) whose Fourier transforms vanish outside [0,1] and satisfy \z(t)\ < m. 



The behavior of A n {m) is such that $ n (m) defined in (241) for < m < 
m n is very close to n over most of this range and increases suddenly as 
m —*■ m n . Consequently, the upper bound B(m) defined in (242) is 
roughly a staircase function as shown in Figure 1. For < m < .62 we 
have B{m) = /3i(m). In Figure 2 a graph of logio \B(m) — 1} is plotted for 
.48 < m < .62. It is seen that only .1% increase in transmission bandwidth 
is required for m < .48, and 10% increase suffices for m < .57. We know 
that /?o(m) > 2 for m > 1, so without Poisson filtering SSBEM is inter- 
esting perhaps only for \z \ < .6. 



X. A POLYNOMIAL PROBLEM 

To each polynomial 



fc-i 



(246) 



we can assign a positive integer v which is the smallest integer >n such 
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-2 

-3 

4 



0.48 



0.50 



0.52 



0.56 



0.60 



0.62 



Fig. 2 — Logarithmic expansion of Fig. 1 for small m. 

that for some choice of a* for k — n + 1, n + 2, . . . , v, the polynomial 
P n ,M) = l+ £ a k { k + ± a k { k (247) 



k=i 



n+l 



is zero free for | f | < 1, The integer v is some complicated function of the 
coefficients a,k, k = 1,2, . . . n. The fact that v is finite is a rather re- 
markable fact that follows from the previous theory. To see this we 
set 



$=re u , r>0. 



(248) 



For sufficiently small r we have P n (re u ) zero-free in the uhp. Then 
taking 

Q(re it ) = log Pnire*) (249) 

= £ b k r k e ikt 
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and 



we have 



Qn(re it ) = Y.b k r k e ikt (250) 

1 



exp\Q n {e it )}=l + 'Lc k e ikt (251) 

i 

where 

c fe = a k , k = 1,2, . . . ,n. (252) 

Now we identify Q n (e it/n ) with z(t) in the previous section where we 
obtained upper bounds on Po(m). We know we can bandlimit exp 
\Qn(e lt )} to obtain a function of the form 

n N 

1 + £ c fc e 1- *' + £ ef*«'** 

1 n+l 

which for sufficiently large N is zero free in the uhp. In particular, we 
have shown that this is possible for 

N < nB(m) (253) 

where B(m) is defined in (242) and 

w-max|Q B (e")|. (254) 

t 

Thus 

v < N < nB{m). (255) 

Given the a,k, or equivalently the bk, for k = 1,2 n, we are inter- 
ested in obtaining a lower bound for v. 
Writing 

*V(f)= II (1-X*f) (256) 

where |X^| < 1 
we have 

logP n ,,(f)= f 6 A r fe , |f|<l (257) 

/e = l 

where 

6k = -7i: (A,)*- (258) 
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Since \\k\ < 1 we have | bk | ^ v/k and hence 

p>max \kb k \, k = l,2,...,n. (259) 

k 

Consider, for example, the case 

P„(fl - 1 + mf», m>0 (260) 

PnAO = 1 + mt n + a n+ i^ +1 + ... a,f. (261) 

We have 

logP».,(fl = mf + 6 n+ ir ,+ 1 + • ■ • (262) 

So for the case (260) we have 

v > nm. (263) 

This inequality is clearly best possible in case m is a positive integer. Now 
suppose m = I + e, where < e < 1/n. Since v is an integer we conclude 
from (263) that v>n + l. However, we can show (see Appendix F) by 
another method that m > 1 in (260) implies u > 2n. Then it is apparent 
from the example 

(with appropriate choice of a 2n ) that v = 2n for 1< m < 2. It is conject- 
ured that this large jump in v at m = 1 also occurs at all integer values 
of m, but we have not been able to show, for example, that m > 2 implies 
j/>3n. 

In order to improve the lower bounds obtained for /3o(m) we are in- 
terested in maximizing the ratio v/n subject to the constraint 



max 

t 



bo + T. b k e ikt 

rfl 



< m. (264) 



For any choice of bk, k — 1,2, . . . n, we are free to choose b so as to 
minimize the maximum modulus of the sum. That is, in the bandwidth 
problem of the previous section we take Q = n and 

z(t) = £ b k e ikt (265) 

k = 

with the constraint (257). Then assuming n < a < n + 1, and v<&<v 
+ 1, we have 



w a , (t) = e b " 1 + Y. a k e iht + £ a k e ikt 

k=l k=n+l 

= e b ° ft (1 -*****), | A* | <1. 



(266) 
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Now we want lower bounds on v implied by (258), subject to (264), giving 
us 0o(m) > v/n. For a given n (large) we would like to choose the bk so 
as to maximize v. For a particular choice of bk we can in principle de- 
termine the minimum value of v required to satisfy (258). In order to do 
this we may assign values to the (n — v) coefficients a*, k m n + 1, . . . , 
v, in (264) and see if it is possible to make \\k\ ^ 1 for k ■ 1,2, ... ,i». 
Recall that the first n coefficients are determined by 



exp 



E b k e ikt 
k=i 



= 1 + one** + . . . + a n e int + ... (267) 



Perhaps a computer study could shed some light on this very difficult 
and challenging problem. 

APPENDIX A 

Proofs of Theorems 1, 2, and 3 

In view of Ex. 2 in Sec. 2.2 it is sufficient to prove Theorem 1 for the 
interval (— °°,0). So we assume first that g(t) is a bounded function whose 
Fourier transform vanishes over (— °°,0) and we wish to show that g(t) 
is the boundary value of a function bounded and analytic in the upper 
half -plane. For this purpose we define 

g u (t) = £° g(s)P u (t - s)ds, u>0 (268) 



where 



We have 



Hence 



Also 



Pu(t)=-—^—- (269) 



\g u {t)\ < J_Js( S )||P u (*- S )|d S . 



sup|£ u (£)| <sup|g(t)|. (270) 

t t 



\\m g u {t) = g{t) for almost all t. (271) 

u— o 

Now all we have to show is that g u (t) is an analytic function of t = t + 
iu. Since the Fourier transform of g(t) vanishes over (— «>,0) we may 
replace P u (t) in (268) by any function of L\ whose Fourier transform 
agrees (for each u > 0) over (0,°°) with that of P u (t). In particular, we 
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may replace P u (£) by the analytic kernel, 

K(t + iu) = — r"F(w)e'" w(t+, ' u >d«, u >0 (272) 

27T «/-» 

where F(u>) = 1 for co > 

= a> + 1 for -1 < co < 

= for u< -1. 

The definition of F(io) for w < is important only to the extent that the 
integral in (272) converges for u > and implies 



f " \K(t + iu)\dt <« for u > 0. 



(273) 



It is sufficient for (273) that e - uu F(a>) belong to L 2 and have a derivative 
inZ/2 [see (5)]. 
Thus we have 

g u (0=— f"g(s)K(t + w-s)ds = G(t + m) (274) 

where G(t) is analytic in the uhp and from (270) and (271) 

|G(r)|£sup|*(t)| (275) 

t 

lim G(t + iu) = g(t) for almost all t. (276) 

u--0 

This proves the first half of the theorem and we may as well write 

G(t + iu) = g{t + iu) = g u (t). (277) 

Now for the second half of the theorem we wish to establish that if g(r) 
is bounded and analytic in the uhp, then 



C~j(t)h(-t)dt = 



(278) 



for all functions hofLi whose Fourier transforms vanish over (0,«), or 
equivalently 



C~j(t)h(t)dt = o 



(279) 



for all functions h of Li whose Fourier transforms vanish over (— «\0). 
To do this we need some lemmas concerning analytic functions belonging 
to Li on lines parallel to the real axis. 

Lemma 1. Ifh(t) belongs to L\ and its Fourier transform vanishes over 

ANALYTIC MODULATION SYSTEMS 539 



(-»,0) then h(t) has an analytic continuation h(t + iu) in the upper 
half -plane u> satisfying 

f" \h(t + iu)\dt< f" \h{t)\dt. (280) 

Proof. We have 

&(*)■ C~ Mt»)e iat da (281) 

with the Fourier integral providing the analytic continuation 

h(t + iu) = f"/i(a;)e l '"<' + '">dco, u > 0. (282) 

Since n(w) = for co < we may write 

h{t + iu) = C fi(u)e- u \ a \e iat du), u>0 (283) 

and hence conclude that 

h(t + iu)= C h(s)P u (t - s)ds (284) 

where P u (t) is the Poisson kernel defined in (269). Then (280) follows 
from (284), since the L\ norm of the Poisson kernel is 1 for every u > 
0. 

Lemma 2. Suppose h(r), r — t + iu, is analytic in the strip <u <b 
and satisfies 



(285) 



Then 



where 



f" \h(t + iu)\dt <°° for < u < b. 
h(t + iu) = — f~ K(sa)et~Q+lu)du t 0<u<b (286) 



fi(a>)= f" h(t)e- iat dt. (287) 

Proof: Defining 

h(u,u) = f°° h{t + iu)e~ iut dt, 0<u<b (288) 

we wish to show that 

K{uyi) - e-» u Ha). (289) 
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We have 

h(t + iu) = lim — f jl - J-^-l fi(u;u)e iut du. (290) 

n-.o, Z7r J-n [ ft 

Also, from the analyticity of h(r), we have 

d Id 

— h(t + iu)=- — h{t + iu), 0<u<b. (291) 

ot i du 

We would like to establish (289) from (291) by differentiating inside the 
integral of (290) but at this point we do not know enough about /i(a>;u) 
to justify the differentiation. Therefore, we will define 

g{t + iu)= C k(s)h(t + iu - s)ds (292) 

where k(t) is a function of L\ whose Fourier transform £(co) does not 
vanish for any argument and such that co£(to) belongs to L\. We would 
also like k'(t) to belong to L\. We may take 

k{t)=-^=e- t2 /*. (293) 

Then g(r) is analytic in the strip and 

f " \g(t + iu)\dt < f " \h(t + iu)\dt. (294) 

We have 

g(£ + iu) = — I ^(co)/i(aj;u)e^ £ da), < u < 6, (295) 
2ir »/-cd 

and since co«(o>) is in Li, 

d 1 /*" 

— g(£ + t'u) = — I i^(u>)fi(o};u)e iwt du), < u < b. (296) 
ot 2ir J-<° 

Now 

£'u (t) = — g(t + iu) = —g(t + iu) (297) 

Ot IOU 

belongs to L\ for < u < b since 

g'u (t) - f /i(s + iuJfe'Cfc - s)ds (298) 

and fe' belongs to Li. Hence the function of t 

d 

— g(t + iu) 
du 
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has a Fourier transform for < u < b. Thus from (295) 

— g(t + iu) = — P"^) — ft(u);u)e iut da, < u < b. (299) 
du 2ir J—*> du 

Since £(a>) ^ for -°° < a; < °°, we conclude from (296), (297), and (299) 
that 

- — ft(co;u) = io>/i(a>;u). (300) 

i du 

Then (289) follows from (300). 

A corollary of Lemma 2 is the following 

Corollary. If g(r), r = t + ill, is analytic in the strip a < u < b and 
satisfies 

f~ \g(t + iu)\dt < *> for a<u<b, (301) 

then 

f " g(t + iu)dt = constant, for a<u <b. (302) 

The corollary follows by applying Lemma 2 to the function g(t + ia + 
U) for arbitrarily small positive «. 

Lemma 3. Ifg(r), t = t + iu, is analytic in the upper half plane u>0 
and satisfies 

f" \g(t + iu)\dt < oo for u>0 

then the asymptotic estimate 

f" \g(t + iu)\dt = 0{e- au \ as u -> « (303) 

implies 

f°° \g(t + iu)\dt <e~ au C" \g(t)\dt for u>0 (304) 

and 

g(£ + iu) = f "gMe^'+^dco for u>0 (305) 

where 

g(o)) = P°° e- iwt g(t)dt = for a> < a. (306) 
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Proof: From Lemma 2 we have the representation (305) so that 

£(u)e-» u = f g(t + iu)e- iat dt. (307) 

Thus 

\§(w)e-™\ < C ~ \g(t + tu)|dt. (308) 

Then (303) and (308) imply 

i(u) = for co < a (309) 

and since ^(co) is continuous, (309) implies 

£(<o) = for <o<a. (310) 

Thus we may write 

g(t + iu) = ?—— C " e-"l^- a l^(a;)e' w d£. (311) 

Then since ° 

1 /* a3 

— I e - u \ u - a \e iat dt = e iat P u (t) (312) 

27r «y— to 

where P u (£) is the Poisson kernel of (269), we have 

g(t + iu) = e~ au f" g(s)e ia ^-«^ P u (t - s)ds (313) 

and (304) follows from (313), and Lemma 3 is proved. 

Now we are prepared to prove the second half of the Paley- Wiener 
Theorem for L„. We have g(r) analytic for u > and 

sup \g(t + iu)\ <M for u > 0. (314) 

t 

Now suppose h(t) is any function of Li whose Fourier transform vanishes 
over (—oo,0). We have from Lemma 1 that h(t) is the boundary value of 
a function h(r) analytic in the upper half-plane u > 0, satisfying 

J \h(t + iu)dt< C \h(t)\dt, u>0. (315) 

Now we consider the function 

f(r) = g(r)h(r) (316) 

which is analytic in the uhp and satisfies 

£~ \f(t + iu)\dt < M f" \h(t)\dt. (317) 
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Then from Lemma 3 

f" f(t)e~ iut dt = for <o < 0. (318) 

In particular 

("°f(t)dt= f°° g(t)h(t)dt = (319) 

which completes the proof of the first (one-sided) Paley- Wiener Theo- 



rem. 



Theorem 2. It is sufficient to prove the two-sided Paley- Wiener Theorem 
for functions g(t) whose Fourier transforms vanish outside [— a, a]. We 
show thatg(£) is the boundary value of an entire function of exponential 
type by defining 

0l(t + iu) = £ ^ ^It^, *. «>0 (320) 

G 2 ( t + m) = S-J^'^ m-^u^ U< °- <321) 
We have 

|Gi(t + iu)\ < e au sup \g(t)\, u > (322) 

t 

\G 2 (t + iu)\< e~ au sup \g(t)\, u < (323) 

t 

and since git) is continuous 

lim G x (t + iu) = g(t) for all t (324) 

u— 

lim G 2 (t + iu) = g(t) for all t. (325) 

u— 

Now we define 

G 3 (t) = f" X(r - s)g(s)ds 
where 

K(r)=^- f"£(«)e i -'rfft> 

27T «/— °° 
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(326) 



and 

#(u>) = 1, -a < a) < a 



= 2(l-£), ,S*S 



2a 



= 0, a > 2a 

K(-o) = £(«). (327) 

Then K(t) is an entire function belonging to Li along each line parallel 
to the real axis. So Gs(t) is an entire function bounded on each line 
parallel to the real axis. Since the Fourier transform of g(t) vanishes 
outside [— a,a] we may replace the convolution kernels in (320) and (321) 
by K(t + iu) since in each case their Fourier transforms agree over 
[-a,a]. 
Thus 

Gi(t + iu) = G 3 (t + iu), u>0 (328) 

G 2 (t + iu) = G 3 (t + iu), u < 0. (329) 

Hence g(t) is the restriction to the real line of an entire function G$(t + 
iu) =g(t + iu) satisfying 

\g(t + iu)\ < e fl l"l sup \g(t)\. (330) 

t 

Now for the second half of the theorem we suppose that g(r) is an entire 
function satisfying 

sup \g(t + iu)\ < e fl l"l sup \g(t)\ (331) 

t t 

and wish to conclude that 

"" g(t)h(t)dt = (332) 



S- 



for all functions h in L\ whose Fourier transforms vanish over (—a,a) 
(and hence over [-a,a]). From the one-sided Paley- Wiener Theorem 
we have 



s: 



g(t)hi(t)dt = (333) 



for all functions hi in L\ whose Fourier transforms are supported on 
(—«,— a) and 



s: 



g(t)h 2 (t)dt = (334) 
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for all function h^ in L\ whose Fourier transforms are supported on (a,°°). 
The difficulty we encounter in establishing (332) is that an arbitrary 
function h in L\ whose Fourier transform vanishes over [—a,a] cannot 
be decomposed as 

h(t) = h-(t) + h + (t) (335) 

where h- and h+ belong to L%, and the Fourier transform of h- is sup- 
ported on (-«,— a), and the Fourier transform of h+ is supported on 
(a,°°). In order to deduce (332) from (333) and (334) we have to ap- 
proximate the test function h in (332) with bandlimited functions hb- 
We may take 



h b (t) - f " bK(bs)h(t - s)ds, 6 > 



(336) 



where 

\& — 



Since 



we have 



sin' 
K(t)=-—^< (337) 

C~ K(t)dt = f " \K(t)\dt = b f " \K(bt)\dt = 1, 

f " \h b (t)\dt < J" |h(t)|rf*. (338) 

Also we may write 

h(t) - h b (t) = f " 6K(6«){/i(t) - fc(« - s)}ds 
which gives 

f " |ft(t) - Mt)|<*t < f " 6K(6s)Mi(s;/i)ds 
where 

m(s;h) = f " |M*) - h(t - s)\ds. 

The function m(s) is called the modulus of continuity of h. It is an even, 
continuous, bounded function of s (see Ref. 3) and 

|w(W-o ^ (342) 

[m(s;h)<2 f" \h(t)\dt = 2\\h\\i. 
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(339) 



(340) 



(341) 



Then writing 

J]_" bK(bs)m(s;h)ds = f " X(Omi (| J h) d* (343) 

J-Vb \b I J\t\>Vb 

it is clear that given c > we can choose b so large (a < b < °° ) that 

f " |fe(t)-M*)|<ft<c. (344) 

Now the Fourier transform of h b is supported on the intervals (—6, —a) 
and (a,6), so hb does have the decomposition 

h b (t) = h-(t) + h+(t) (345) 

desired in (335). This follows from the existence* of a function K aib {t) 
in Li whose Fourier transform satisfies 

R a ,b(w) = 1, a < a> < b 

= 0, a; < -a (346) 



so that 



and 



We have 



h+(t) = f " h b (s)K a>b (t - s)ds (347) 

llfc+lli< ll/ifclli-l|K ai6 |i. (348) 

h-(t) = h b (t) - h+(t). (349) 



So h_ also belongs to Li. 
Returning to (332) we have 

j~j{t)h{t)dt= fjj(t)h b (t)dt + C~j(t)[h(t) - h b (t))dt. 

(350) 
Since 

£~ g(t)h b (t)dt = £" g(t)h-(t)dt + C " g(t)h + (t)dt (351) 

= + 



* See Appendix B for a good choice of K^U). 
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we have 

I f " g(t)h(t)dt < sup \g(t)\ J]_" |M«) - h b (t)\dt 

<€sup\g(t)\. (352) 

Since we may choose b sufficiently large to make e arbitrarily small we 
conclude that 

f" g(t)h(t)dt = (353) 

for all h in Li whose Fourier transforms vanish over (-a,a). This com- 
pletes the proof of the two-sided Paley- Wiener Theorem. 

Theorem 3. Here we wish to show that if g(r) is analytic in the uhp and 
bounded on the real line as well as every line parallel to the real axis in 
the uhp, then the asymptotic estimate 

sup \g(t + iu)\ = 0{e- au \ as u -* «> (354) 

t 

implies 

sup \g(t + iu)\ < e~ au sup \g(t)\, u > 0. (355) 

t t 

Now if x is any real number and y is any positive number, the func- 
tion 

. , , y \g(r)e- iaT - g(x + iy)e-'°(*+'»} 

h{T) = (<35b) 

7T (r-x) 2 + y 2 

where we think of x and y fixed, is analytic in the upper half-plane u > 
and satisfies 

C" \h(t + iu)\dt<«>, u>0 (357) 

and 

C~ \h(t + iu)\dt = 0(l) asu^o. (358) 

It follows from Lemma 3 that 

f" h(t)dt = 0. (359) 
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Hence 

g(t)e~ iat dt 



y /**" s(t)e~ l dt 

77 7Tr- 2 = S(x + iy)e- ia{x+iy) - (360) 

Thus 

|g(z + iy)| < e~°y sup |g(t)|. (361) 

t 

This inequality holds for any -°° < x < °° and any y > 0. Then (355) 
follows from (361), if we note that (355) holds trivially for u = 0. 

APPENDIX B 

Reproducing Kernels of Small L\ Norm 

We would like to find a kernel K a ,p(t) of minimum L\ norm whose 
Fourier transform satisfies 

R-a.&W) - 1, < 0) < a 

= 0, a> > (362) 

where < a < ft. 

Replacing K at p(t) by l/aK a ^(t/a) we see that the minimum norm 
is a function of /3/a, or if we like, a function of a/(/3 - a). It is sufficient 
to consider functions K\(t) whose Fourier transforms satisfy 

# A (a>) = 0, « < (363) 

= 1, 1<«<1 + X 

where we make the identification 

X = T^- (364) 

p — a 

We will not treat the minimization problem here. Instead we give a 
construction for a particular function R\(u)) which can be shownt to be 

the solution for the case X = n,n = 1,2,3, The construction provides 

an interpolation between the minimal norm values in case n < X < n + 
1. 

We write 

X = n + 6 (365) 

where n = [X] is the largest integer contained in X and < 6 < 1. Then 
we set 

K x (t) = 2ir\F x {tW (366) 



f The details will be given in a future paper. 
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where the Fourier transform of F x (t) vanishes for negative argument, 
and otherwise is defined by 

P x (ta) = a k , k <w<k+l <n + l 
= a n+h n + l<a><n + l + 
= 0, a>n + l + d = \ + l (367) 

where the a* are defined by 

- 7 L= = £a k z k (368) 

Vl-2 

ak = ~kT~ fci 



(H 



r(i)ra + « 

We then have 

^M= JJPxWhia-xWx 
which is a piecewise linear function satisfying 



(369) 



(370) 



m-l 



&x(m) = L a fe a m _ 1 _ fe = 1 (371) 

k=0 

for m = 1,2, . . . ,n + 1. Thus 

£ X M = 1 for l<co<n + l 

= a> for 0<co<l. (372) 

For rc + l<co<n + l + 0the convolution in (370) is independent of the 
definition of P\(x) for x > n + 1 + 0; i.e. 

# x (w) = tfn+iM for O<a><n + l + 0. (373) 

So 

# x (co) = 1, n + l<(c<n + l + 6. (374) 

Thus 

# x (a>) = 1 for l<a><n + l + d 

= a; for < a> < 1. (375) 
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We have 

M(X) = f " \K x (t)\dt = 2vr §2 \ F M\* dt 

= f~ |Fx(")| 2 da> 

= IoI + ^ +1 (376) 

o 

which is a piecewise linear function of X, i.e., 

n(X) = (n + 1 - XMrc) + (X - n)M(n +1), n < X < n + 1. (377) 

We have 

5 89 

M (0) = l, m(D=t, m(2) = -- (378) 

4 64 

The remainder of this appendix is devoted to estimating n{n) for large 
n. We need an upper bound. 
For convenience we set 

^ ) = Y^) (379) 



M(n)=- L 7 2 (fe)- (380) 

First we estimate y(x). We have the representation for the Beta 
function 



T(x)T(y) 



= C ty-Hl -tY-hlt (Rex >0,Rey > 0). (381) 



T(x + y) 
Then 

7 (x)=— f * -L. dt. (382) 

vtt Jo Vt(l - t) 

Setting 1 — t = e~ s we obtain 

1 « mp -s(x+l/2) 

1 ~a> p -s(x+l/4) 

= -7= I 7 ; da. (383) 

ViJo V2 sinh s/2 
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Since 

l-e~ s <s 
and 2 sinh s/2 > s we obtain from (383) 



J <7(*)<— 4=' 084) 



y/, + \ Vx + i 



2 4 

We can obtain a simple upper bound for n{n) from (384). We have 

n(n) = l + - ± 7 2 (k) 

7Tfc=l 
1 n 1 

<1+-E • (385) 

4 
Since £ -1 is convex, we have 

J*T+ldt 1 
— > , T > (386) 

r t 1 



T + - 
2 



and thus 



1 /»n+3/4/Vf 1 / 4 n \ 

M(n)<l + - f f =l + 1 log(l+^Y <n>l). 

7T «/3/4 t 7T \ 3 / 



(387) 



Since log (1 + 4X/3) is a concave function of X and since n(\) is piecewise 
linear between integers, we conclude from (387) that 



/z(A)<l+^log(l+y), A>0. 



(388) 



A sharper estimate of m(X) for large X is obtained as follows. We are 
interested in the constant term in the asymptotic expansion. 
From (383) we have 



\ 2/ v^ Jo Vl-e~ s 
Then from the convolution theorem for Laplace transforms, 

y 2 (x-~) = f e- sx <p(s)ds 

where 



7T «/o vl — e~ l 



dt 
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(389) 



(390) 



(391) 



Setting e ( = u in (391) we obtain 

M = l f ' ' , ^ (392) 

TT Jv Vl — WVU—VVU 

where 

u = e~ s . 
Then with the substitution 

l-u= (l-v)t 
(392) becomes 

« {s) = ; Jo ^PiVRr^' (393) 

Identifying <p(s) with the hypergeometric function which has the 
representation, 

2 F l (a,b\c;z) = J ic) M C t b ~Hl - t)«-*-i(l - **)-«(** 

r(6)r(c - o) jo 

{Rec>Reb> 0) (394) 
we see that 



'1 1 

Thus 



ipto-«Fi(|,|; !;!-«-•)■ (395) 



r°(* +i ) 



(396) 



and 



In 1 /»•» 1 — p-*(n+l) 

m(») — Z 7 2 (*0 = " ( — ; — e--/M*)<fe 

7rA=o t Jo 1 — e s 

= I r"{ € -*/2_ c -5(n+3/2)} g^J _ 1| 

7r Jo 1 — e~ s s\ 

7T JO 



e -s/2 _ e -s(n+3/2) 



= ilog(2 n + 3) + i f" e -^.|-2<!>— IU 

T 7T Jo [I — e~ s s\ 



--s 

7T JO 



e -s(n+3/2) 



' **--±i* (397) 



1 — e s s 
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From (395) and the series 



2 Fi(a,b;c;z) = £ 



ia)k(b) k z k 



it is clear that 



and hence that 



k=o (c) k k\ 



(f)(8) > <p(0) = 1, s > 



(398) 



(399) 





1 


*> 1 > „ 

- e s s 








(400) 


Thus from (397) we have 






/i(n)<-log(2n + 3)+- f°V« /2 

IT IT JO 


<p(s) 


ll 

s 


ds 


(401) 


1 -e _s 


and in fact 






lim 

n— ►<*> 


M (n)--log(2n+3) 1 


= - f°V«/2 

7T JO 


[ v(s) 
1 -e _s 


_1 

s 


ds. 


(402) 



Now we can evaluate the integral in (402) by an indirect route. 
We have 



1 5 1 n 

v(n)—Z— 7 = E 

■K K + 1 



(1/2)* (1/2)* 



-if: 

7T 



2fM _ 



y 2 (k) 



k!k! v(k + l). 

1 l 



fe + 1 



(403) 



From the estimate (384) we see that the sum on the right converges 
as n -* «; i.e. 



1 n i 

lim M(n)--Ir-T 

7T « + 1 



1 " f 1 

= -Z T 2 (*)-r— : 

X [ « + 1 



(404) 



Now we can write 



7T 



yHk) __i_| . lim [ £ flawk»k,» _ 



fe + U «-i U-b fete! 



= lim 

X — 1 



1 1 



>MH :1; *)-> g r^) (406) 
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Then using (394) we have 

1 



Vt(l -0(1 -xt) 1 -xt 



irA = lim ( 
x— l Jo 

= lim p_*_L'IZiL_ 1 | 

,-i Jo l-xt\Vt(.l-t) 

-lim f'-^-fi-a + ita f — *_ _ 

x-iJo 1-xtlVt J x -*i Jo (l-xt)Vt VI- 1 



Vl -xt 



-1 



■X' 



dt 



(l + VI) Vt 



+ 



lim f 1 - 

x-lJo (1- 



dt 



xt)Vt 



Vl -xt 



- 1 



= 21og2 + 21og2. (406) 



Some care is required in evaluating the last limit. An alternative way 
of obtaining the result is worth noting, as it makes use of an interesting 
series obtained from (396). A change of variables gives 



r2 (* + 2) 

(1/2M1/2)* 



dt 



-i 



o I o 



k\k\ 



t k \dt 



x *(x + l) 



ll 



+ 



2! 


2 2 
2 


x(x+ 1)(jc + 2) 



+ 



(407) 



Then setting 



( (1/2)/ 
1 k\ 



^ ( X + I) 



r (|) r<* + 1) 



2 l 
= -F(fe) 

7T 



we have 



F(x-l)=-t-TV—F(k) (408) 

\ 2/ 7Tfe= (x)/ ?+1 

which is an interesting formula. In particular (F(0),F(1),F(2), . . . ) is 
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an eigenvector of a certain infinite matrix. Also an interesting series for 
ir is obtained by setting x = n + 1/2 (large n) in (408). Returning to (407) 
and recalling a* = (l/2)k/k\ we have 



xrHx) r 2 (x + i) af | 



2! 



r2 H) * r! (* + i) 



x + 1 (x + l)(x + 2) 



(409) 



Then using the series 

-- C tt-t) x ~ l dt= C 1 (l t -t) x (l + t + t 2 + ...)dt 

x Jo Jo 



+ 



+ 



2! 



x + 1 {x + l)(x + 2) (x+l)(x +2)(x + 3) 
we may write 



+ . . . (x > 0) 



1 I T 2 (x + D T 2 (l) 

*V(-i)"'-© 



Since 



= 1 + 



x + 1 



2 H 

7T 



+ 



2! J 9 1 

(jc + l)(x + 2) T 2 2tt 



+ ... 



fe! 



(x + l) k 



al- 



+ ... (410) 



»(* + i) ,(* + i) 

we may let x -»■ with the result 



t . - r , il d r 2 (x + 1) 

fe=i l «ttJ dx 



r 2 



K) 



= -log2 (411) 



x=0 



and since ao = 1, this sum is the same as the sum on the right in (404). 
Hence the limit in (404) is 



lim 



1 " 1 
M(n)--L 



= -log2 

IT 



TTok + 1 
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(412) 



and since 
lim 



« 1 

£ - - log n 
1 n 



= C = 0.577215. . . (Euler's constant) (413) 



we have from (412), (413), and (402), 



lim 



M(n)--log(2n + 3) 

IT 



■K JO 



,-s/2 



1 — e -s s 



ds 



= ^log2 + ^ 

■jr x 



Then from (401) and (414) we have 



1 3 C 

n(n) < - log (2n + 3) + - log 2 + - 

7T 7T TV 



and by the same argument used in establishing (388), 

M(X) < - log (2X + 3) + - log 2 + - , X > 0. 

IT TT IT 

We find from (397) and (414) that 

irn(n) ~ log (2n + 3) + 3 log 2 + C 

3 43 7 



(414) 



(415) 



(416) 



2(2n + 3) 48(2rc + 3) 2 16(2rc + 3) 3 

+ 0(n" 4 ). (417) 



For comparing the estimates (389) and (415) we have 



1 + - log - = 1.0915720476 
■k 3 



- {4 log 2 + C\ = 1.066275853 



(418) 
(419) 



and for use in (415) and (417) 



1 .« , n ^ 2.656657207 
- 13 log 2 + C( = 

IT 7T 

= .8456402533 (420) 
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In the following tabulations the estimates (388) and (415) and the as- 
ymptotic formula (417) are compared with the true value of n(n). 



n 


m(m) 


Asymptotic formula 


Error 


1 


1.25 


1.249927097 


7.29 (-5) 


2 


1.390625 


1.390607982 


1.70 (-5) 


3 


1.48828125 


1.488275479 


5.77 (-6) 


4 


1.563049316 


1.563046870 


2.45 (-6) 


5 


1.623611450 


1.623610248 


1.20 (-6) 


6 


1.674500465 


1.674499809 


6.56 (-7) 


7 


1.718379259 


1.718378872 


3.87 (-7) 


8 


1.756944605 


1.756944363 


2.42 (-7) 


9 


1.791343941 


1.791343782 


1.59 (-7) 


10 


1.822389342 


1.822389234 


1.08 (-7) 


n 


Upper bound (388) 


Upper bound (415) 




1 


' 1.269703286 


1.357940252 




2 


1.413574619 


1.465042692 




3 


1.512299999 


1.545038559 




4 


1.587544884 


1.608914025 




5 


1.648359655 


1.662088992 




6 


1.699398305 


1.707639405 




7 


1.743372924 


1.747480071 




8 


1.782003284 


1.782884290 




9 


1.816448738 


1.814741844 




10 


1.847528028 


1.843699061 




APPENDIX C 






Estimates for M„ 







M n = 



E cot 



kir 



n+1 i< k <n 2(n + 1) 

k odd 



(421) 



In order to express the sum as an integral we note first that 



■K COt 7TX = — log Sin 1TX 

dx 

d 

= — -log 



dx T(x)T(l-x) 
= t(l-x)-yp(x) 



where 



nx) r(x) Jo I logt i-t 



dt, x>0. 
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(422) 



(423) 



So 



7T COt TTX 



=j>-— ^ 



(424) 



■X 



U-u(x-l) _ e uxj 



du, < x < 1. 



Then for < 6 < 1/p, where v is an odd integer, 

( e -u0 _ e -(H-2)u0) 



-2u0 



7T0 £ COtfe7T0 = j e" _, 



k odd 



1 _ e 2u0 



1 -e~ u 



du (425) 



or 



■kO £ cotkitd 
k=\ 

k odd 



-I 



(e- t -e- ( " +2)f ) (e-'-e*) ) 



,-at 



dt 



1 -e" 2 < l-e" 2£ J 1 -e~ at 

^ <-.-»* r (l- e - W ) t )2 e <-o), 

Jo 1-e" 2 ' Jo l-e" 2£ l-e- a < 

where a = 1/0. For n odd we take v = n,a =2(n + 1). Then 



,-* _ e -i.n+2)t 



M _ e -(M+l)tj2 e -(n+2)t 



M " = ;Jo l- e -» dt ~*Jo l-e~* l-e- 



,-t _ e -(n+2)t 



4 -,a> e -« _ e -(n+in 2 /» m l - e L at 

= * Jo l-e" 2 ' f ~7r(n+l)Jo . i t l + e £ 



2(n+l)( 

l_^e dt 

sinh 



dt 



n+l 

n odd (427) 

The first integral is just the sum of the reciprocals of the odd integers 
from l through n. We have 

e -t _ e -(n+2)t 



I 



t 



dt = log (n + 2). 



So 



-t _ e -{n+2)t 



J»°° e~ L — e v 
o l - e" 



2< 



dt = log (n + 2) + 



r«- 



l-e- 2 ' t 



(428) 



-±U 



-(n+2)t 



-X Vl " + 



-- dt. 



(429) 



l-e" 2 ' t 
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We have 



Jo ' il-e-2' t. Jo il-u 2 logu 



where 



Thus 



log 
du 

du 



du 



. f'LJ_ + _!_ + _L_ 

Jo [ 1 — u 1 + u log u 

»/o 1 + w Jo 1 1 — u log u 
= log 2-^(1) 

= log 2 + 7 



7 = .577215. . . (Euler's constant). 



(430) 
(431) 



M n = -log (n + 2) + - (log 2 + 7 ) 

7T 7T 



7T Jo 



d* 



^-T f 

7r(n + 1) Jo 



l-e~ 2t t 
l-e- 1 dt 
t l + e« 



sinh 



n+ 1 



We may write 



I 



e -(n+2)t 



1 - e~ 2t t 



(n odd) (432) 



dt 



= r " e-^+D' [ — - — - - — i 

Jo Isinht t 

= f" e -(«+i><(— A M«ft + f" e -(»+l)«(l-e-*) — 

Jo Isinhi £j Jo t 



-4-r.- _j — 

n + 1 Jo 1 t 

I sinh 



«+'l* + toi a±« (433 ) 

M + 1 



n + 1 



Also 



1 p~ l-e- f 
n + 1 Jo 1 + e l 



dt 



sinh 



t 



n+ 1 
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n+ 1 



"n + lJo 1 + e' t t 

I sinh 



dt 



n + 1 



l-e-'dt 



; — (434) 

o 1 + e' t 



and (Ref. 9, p. 327) 



l-e-'dt 



J»«° l -e ' dt 7T 

o l + e £ t " ° g 2" 



(435) 



Thus we find from (432)-(435), 

2 2 / 4 \ 

M„--log(n + l)+-(log- + 7) 

+ 4 f" 

7r(n + 1) Jo 



n+1 



dt 



sinh 



t I 1 + e 1 



n+1' 



n odd. (436) 

Clearly the integral in (436) is positive. In order to obtain the asymptotic 
series for M n , n odd, we use the generating function for the Bernoulli 
polynomials (Ref. 1, formula 23.1.1) 



to obtain 



fpXt <*» f-k 

e l - 1 k=o k\ 



— - — £ * (I) S w<*r 

t ^=o \2/ fe! 
2 sinh - 
2 



(437) 



(438) 



where (Ref. 1, formula 23.1.21) 

Bk (2) = " (1 " 2l ~ k)Bk ' B * sB *< 0) " 



Thus 



e 



t* 



_ J*-2)B*- 
sinh t a=o «' 



= - £ (2 fe -2)B*^ |t| <T. 



(439) 



(440) 



We have (Ref. 1, formula 23.1.19) 

B 2k +i = 0, k = l,2,... 



(441) 
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and (Ref. 1, Table 23.2) 



So 



Bo-1, S i-"2' fi 2--, fi 4-- 30 

S6 "42' B8_ 30' Sl0 "66' Bl2 " 2730' 

Bu = I • (442) 

6 



7-^-=E(2 2 "-2)^^ (443) 

t sinn £ k=i (2ft)! 

= "-— t3 + ^-t 5 + ... |i|<7T. 

6 360 15120 
From Ref. 9, p. 325, we have 

J, * r 2A-l -.2A 

-*— dx=(2 2 *-2)|B 2fe |-— , k = l,2,... 
o e* + 1 4ft 

Also, (Ref. 1, formula 23.1.18) 

B 2 „ = B 2 „(0) -(-!).« ^tZ^, -1A... (445) 

SO 

fi 2fe =(-l) A+1 |fi2fe|. (446) 

Thus we obtain the asymptotic series 



— f " 

n + 1 Jo 



n + 1 1 j dt 



t ( 1 + e' 
sinh 



- (-l) fe+1 (2 2 * - 2) 2 Bl k / 7T ^ 



(444) 



n+ 1 

- (-l)fe+l(2^ - gjjggfc / 7T y 

fc Ti 4ft(2ft)! \n + 1/ 

72 Vn + 1/ 43200 Vrc + 1/ 
and hence 

2. . _ 2 /. 4 \ . 7r 49 7T 3 



Z Z / 4 \ 7T 4i7 7T 

^n— log(n + l) + -(log ; + 7 j + __-_ io8oo(w + ^ 
+ ...+ 



(-!)»+. (2»-2)Ug * /_J_) M + ... (nodd) . (448) 

7T ft(2ft)! Vrc + 1/ 
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Now we would like to show that 



M„>-log(n + l) + -Aog- + Y) 

7T 7T \ 7T / 



<-log(n + l)+-(log- + 7 )+— — A 
ir ir \ ir / 18(n H 



18(n + l) 2 



>-log(n + l) + -(log- + 7 )+^-^ 
7T T \ IT I 18(rc H 



49tt 3 



+ l) 2 10800(n + l) 4 
etc., (n odd) (449) 

i.e., that the error in truncating the asymptotic series has the same sign 
as the next term of the series. To do this, we show that for t > 



>0 



t sinh t 

<*- 

6 

6 360 

<£_JL ,, + _«_,. 

6 360 15120 
etc. 



We have (Ref. 9, p. 23) 

t 



= 1 + 2t 2 V (-l) k 

sinht fc t*i ' t 2 + k 2 ir 2 



or 



t 2 sinh t 



t l co 1 

= 2Z(-ip „;. 2 2 . 



Now consider the polynomial 



P 2 n(t;k) = 2 



\ik-K/ 



I \2n+2 



t 2 + k 2 TT 2 ' 



We have 



(450) 



(451) 



(452) 



(453) 



P2n(t)= L (-l) k+1 P 2n (t;k) 



1 

t 2 



1- 



sinh t 



_ Of — -Wn+lt 2n+2 



2(-l) n+1 t 



L 



-■wk+i 



(-D 



(454) 



k=i(kir) 2n+2 (t 2 + k 2 Tr 2 ) 
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or 

1 

t 2 



^l 1 -^}-^) 



= 2(-l)" +1 £ 2 "+ 2 £ 



-1\A+1 



1) 



(455) 



k=l(kTr) 2n+2 (t 2 + k 2 TT 2 )' 

It follows that P2n(t) is a polynomial of degree 2n which agrees with the 
first (2n + 1) terms of the Taylor series of t~ 2 \l — t/sinh t\. The sign of 
the difference in (455) is (— l) n+1 , as the sum is clearly positive. Then 
(450), and hence (449), follows from (455). 

For n an even integer, the asymptotic expansion of M n is not so readily 
obtained. In this case we set v = (n — 1) in (426) and keep a = 2(n + 1). 
Thus 



Mn = ~ 



■K JO 



,-t _ e -(n+l)t 

l-e~ 2t 



dt 



4 r™ (1 - e~ nt ) 2 e~ (n 
"ttJo l-e~ 2t l-e~ 



e -(n+3)t 



2{n+l)t 



dt, n even. (456) 



We have 

i a, Q—t _ g— (n+l)t 



I 



l- e -2t 



+ 



dt = log (n + 1) 
r°°{c- f -e-<" +1 ) f } 



l- e -2t 



dt. (457) 



Now we would like to express the second integral in (456) as an asymp- 
totic series in (rc + l) -1 . For convenience we set e~ l = x. Then 



(1 



_ a -nt\2 



)2 (l- X n)2 { 1 _ x a+l_ x n( 1 _ x )}i 



I _ e -2(n+l)t i _ x 2n+2 (1 _ x n+l) (1 + x „+l) 



= 1- 



x n (l -x)} 



1 -x n+1 



l-x n+1 x n (l-x)\ 



l + x" +1 



l-x n+] x n (l~x) x n (l-x) X 



l + x n+1 

2n(l _ x )2 



1+X n+1 1 + X n+1 1+X n+1 l-x 2n+2 

Therefore 

(l_ e -rat)2 e -(n+3)t 



I 



1 - e~ 2t 



dt 



J*°° 1 
o 1 



I _ e -{n+l)t e -{n+3)t 



I - e -2(n+l)t~ J Q 1 + e -(n+l)t I _ c -2« 

e~ nt {l — e -t ) e~( n+3)t 



dt 



-2 r e - 

Jo 1 

+ x 



+ e -(n+l)t l — e -2t 

e -2nt(l _ e -t)2 e -{n+3)t 

I _ e -2(n+l)t I _ e -2t 



dt 



dt 
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I 



I - e -(n+l)t fa 



o l + g<» +1 >« e 2t - 1 

J, co c -2(w+l)t <jt 

o l + e-< n+1 ) f e'+ 1 

1 - e -2(„ +1 , tl + e - t dt - <««») 

We can combine a part of the second integral with the last by noting 
that 

= tanh 



e' + l 2 2 2 

Then 



e -2(n+l)t ^£ ^» e -3(n+l)t £ 

tanh - dt 



J»» e~' WITtw dt /"» c~ ai " 

o \ + e - (n+l)t e'+ 1 Jo 1-e" 5 

= -I :<*£ + -; : -tanh— -dt. 

(n + l)Jo 1 + e-' (n+1 Jo 1-e" 2 ' 2(n + 1 



By making use of (430) and (435) we obtain 
M n = -login + 1) + - (log- + y) 

■K IT \ IT / 

ir Jo \l-e- 2t t 

l-e-' \n + l 



2 /»- 1 - e l \n-t 

7r(n + l) Jo 1 + e' ) t 



rft 



2t 
exp 



7r(n + 1) Jo 1 



,-2t 

dt 



4 r~ e" 2f , t 



4 /* ea e~ t 

— I — tanh dt, n even. (459) 

7r(n + l)Jo l-e~ 2t 2(n + 1) 

From (437), (441), and (442) we have 

2 1 Bht k ~ l ■" Bobt 2k ~ 1 



e 2t - 1 t fc ti fe! fetTi (2fc)! ' 

(460) 

r^-r^l.^w 1 ' i*i<- (461) 
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i.e., 

Thus 
2 



2 -Ll a ^_-i + l=E2-^. (462) 



l- e -2t t e 2t -l t k =i (2k)\ 



l-e-t \n+l 



w(n + 1) 



J»°° 1 -e~ 
o l + e f 



t / 2i 

exp 



dt 



\n+l/ 



-- f%-<»+i>t( ^:-- di 

7T »/0 



l-e~ 2t t 



2 r~ l-e- 

" tt(m + 1) Jo 1 + e< 



ir(n + 1) »/o 



l-e- 1 \n+l 



exp f } - 1 



-1\ dt 



\n+l/ 



M + l 



, 2 /•- fl_-e-' 

tt(m + 1) Jo 1 1 



-lldt 



7r(n + 1) Jo 1 1 + e< 



+ e' 



\n+ 1 



-c-' dt 



* / 2t 

exp 



-1 \dt 



\n + 1/ 
4 /»* e -t 

7r(n + 1) Jo 1 + e< 



dt. 



Then from the above and (459) we have 

2 2 / 4 \ 

M n = -log (n + 1) + - ( log- + t) 



4 p°° [ n + 1 2 

ir(n + 1) Jo j t i it \ 

( exp Cm) - 

7r(n + 1) Jo 1 - e- f 4(n + 1) 



-1 



dt 



1 + e' 



-t t 

tanh — — —7 dt, n even (463) 
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which is a form suitable for asymptotic expansion. We have 

v2fe-l 

tanhx= L 2*H2 2k ~ DB 2k — — 
fe-i (2ft)! 



(464) 



t k -= I t k (e- t -e- 2t + e- 3t + ...)dt 

o l + e< Jo 

= k\(l-j^ttk + D (465) 



where 



Also 



k=\k n 



(466) 



f °° — -t k dt = f " i*(e-' + e~ 2t + e~ 3t . . .) dt 

Jo 1 - e~ l Jo 



From (460) and (465) we have 



4 /•- I n + 1 

?r(n + 1) Jo j t 



exp 



(— )- 



4 - J3 2fe 22fe(2fe - 1)! 
t*-i (2ft)!(n + l) 2fc 



-1 



= ft!f(ft + l). (467) 



dt 



1 + e' 



(i-^t)««) 



and from (467) and (464) 
7r(n + 1) Jo 1-e"' 



tanh 



-l£ (22 *_ 2)^-^j (468) 

7rfe=i ft (ft + l) 2 * 



dt 



4(n + 1) 
B 2 k 



-- T 2 2fe (2 2fe - 1) 

*k-i (2k)\A 2k -Hn + l) 2k 



(2k - 1)! f(2ft) 



7T/e = l ft (n + I)** 
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From (445) we have 



B2*=(-l)* +1 2-^-f(2fe). (470) 



Adding (468) and (469), using (470), we obtain 

2 2 / 4 \ 

M n log (n + 1) + - (log- + 7) 

7 7T 127 7T 3 



36 (n + l) 2 21600 (n + l) 4 

+ ; ( - 1)4(244 - 2) Ai^ + - ™- (471) 

In the same manner as before, we can establish that 
M n <-log(n + l)+-/log- + 7) 

>-log(rt + l)+-(log- + 7)-^— ^ 
■k 7r \ 7r / 36 in + 



<-log(M + i) + -(log- + 7)-^— fTTi + 
7r 7r \ 7r / 36 m + l) z 



(a + 1) 2 

127 7T 3 



(rc + l) 2 21600 (n + l) 4 
etc., for n even. (472) 

In this case we have two functions to consider in the polynomial ap- 
proximation problem. First we note that 

^~ l + 1 = cothx (473) 

and 

j W r co 2 

coth x = — log sinh * = — log x n (1 + TTTM 
dx dx b k =i \ k 2 * 2 / J 



Clearly, for x > 0, 



= 7 +2 * s 2 rti 2 • < 474 > 

x jt=ix^ + ft 2 ir 2 



^T-J«> (475) 



and 



--"-•I l = cothx--=2x E — 7TT >0. (476) 



e 2x — 1 X X fc = lX 2 + ft 2 7T 2 
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Now we have 








coth x 


1 








X 9 f 


1 




X 


k = \X A 


+ fc 2 7T 2 




Defining as before 








1 - 

P2n(x;k) - 


/_X_\2n+2 
\ikir/ 

t 2 +fe 2 7T 2 


fe = 


1,2,3, • • • 



(477) 



(478) 



and then 



Q 2 n(*) = 2 LP 2 „(*;fc) 



coth x 

X 



2n+2 



X 

we have 



coth x 

x 



-^ E 

7r/ /e = l/2 2 " +2 (jC 2 + fe 2 7T 2 ) 



(479) 



2n+2 



(480) 

So Q2n(jc) is a polynomial of degree 2n which agrees with the first (2n 
+ 1) terms of the Taylor series of x _1 {coth x — x~ l \ and the sign of the 
difference is (— 1)" +1 . 

For the second function we have 

tanh x = 4~ log cosh x = 2x E . (481) 

4 
Now we define 

2x \2n+2 



1 



\i(2k - !)*■/ 



P2n(x;k) = ^ '-^-j-, k = 1,2, . . . (482) 

x 2 + (2fe-l) 2 v 
4 
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and 



qtn(x) = 2 £ P2n(x;k) 



tanhx 



(9y\ 2ra+2 » 
-) E 
7T/ fe = l 



( 2 fe - l)2n+2 | x 2 + (2fe-l)2 — 



(483) 



Then the alternating sign of the error 

tanhx 



follows; i.e., 



-<72n(*) 



tanh x < x 

. x* 
> x 

3 

r 3 o 

<x-- + —x 5 
3 15 

etc., for x > 0. 



(484) 



The inequalities (472) then follow from (463), (476), (480), and 
(484). 
Then for n even or odd we certainly have 



M n > - 



;>-iog<» + i) + -(iog- + 7) 

7T 7T \ IT / 



36 (n + l) 2 



or, giving away a little, 

2 



z z / 4 \ 

M n > - log n + - ( log - + y ) , n > 1. 

7T IT \ 7T / 



(485) 



(486) 



Obviously (486) is true for n odd since for n odd we may replace n in 
(486) by (n + 1), and for n even we consider 



log(« + 1) -^-^— = logn + log (l + -) -^-^ 2 
72 in + lr V n/ 72 m + 



(n + 1) 



2 ' 



Now 



> I {l-t)dt = x , x>0. 

o 1 + 1 Jo 2 
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So 



log / H i\.l^ > i_i.i >0 for „> 2 . 
e \ nl 72 (n + l) 2 n 2rc 2 rc 2 

Thus (486) is valid for n > 1. 

We note that the actual computation of M n is simplified, particularly 
for n odd, by making use of the identity 



cot0 = 



1 + cos 20 
sin 20 



Then for = w/2(n + 1), 
2 



M n = 



M + 1 l<k<n 
k odd 



£ cot kO = 



1 + cos 2kd 



n + 1 i<A<n sin 2&0 

ft odd 



For n odd, fc = 1,3,5, . . . ,n, 



So 



2/ztt 2(n + 1 - k)ir 

COS — ~ ~ 7 = —COS 



2(n + 1) 



2(n + 1) 



2/ztt . 2(n + l-/z)7r 

sin —. — = sin 



2(n + 1) 



2(n + 1) 



(487) 



We find 



M n = 



fl + 1 l<fe<n . &7T 

ft odd sin 



n odd. 



n+ 1 



M 2 = 



M x = l 

2\/3 



= 1.1547005. . . 



\/3 3 
M 3 = V2 = 1.4142135. . . 



(488) 



M 4 = - 



1 2 

+ 



. TT . 27T 

sin — sin — 
5 5 J 



R • 27r 
5 sin — 

5 



[1 + COS 7r/5] 



= 1.5216904. . . 



M* = - = 1.6666666. . . 
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For use in the asymptotic formulas we have 
- = 0.636619772367 . . . 

log - = 0.241564475270 . . . 

7T 

7 = 0.577215664901 . . . 
log - + t = 0.818780140172 . 

■K 

- /log - + 7) = 0.521251626 

7T \ IT I 



18 ~ 


0.174532925 . . . 


7_7T_ 

36 


0.610865238 . . . 


49tt 3 
10800 


= 0.140676625 . 


1277T 3 

21600 


= 0.182305423 . 


The inequalities (449) and (472) give 



Mi > 0.96252282 . . . 

< 1.00615605 . . . 

> 0.99736376 . . . 

M 2 < 1.2206499 . . . 

> 1.1527760 . . . 

< 1.15502669 . . . 

M 3 > 1.40379402 . . . 

< 1.41470233 . . . 

> 1.41415281 . . . 

M 4 < 1.54585162 . . . 

> 1.52141701 . . . 
< 1.52170870 

M 5 > 1.66192113 . . . 

< 1.66676926 . . . 

> 1.66666071 . . . 
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APPENDIX D 
Proof of Theorem 5 

We are given 

w(t) = e 2 «> (489) 

where the Fourier transform of z(t) vanishes outside [0,Q]. Also the 
Fourier transform of w(t) vanishes over (a,b) where < a < b and 

b - a > fi (490) 

and we wish to show that w(t) = constant. 
We may write 

w(t) = g(t) + h(t) (491) 

where the Fourier transform of g(t) vanishes outside [0,a] and the 
Fourier transform of h(t) vanishes over (— «>,6). We have 

w'(t) = z'(t)e*M = z'(t)\g(t) + h(t)} (492) 

= g'(t) + h'(t). 

Now the Fourier transform of z'(t) vanishes outside [0,12] and the Fourier 
transform of g'(t) vanishes outside [0,a]. By Corollary 2 of Theorem 2 
the Fourier transform of z'(t)g(t) vanishes outside [0,a + Q], The Fourier 
transform of h'(t) vanishes over (— °°,b) and by Corollary 1 of Theorem 
1 the Fourier transform of z'(t)h(t) also vanishes over (— °°,b). Thus if 
K a b is any kernel of L\ satisfying 

C* K a ,b(t)e- iut dt= 1, 0<co<a (493) 

= 0, w > b 



we have 



f" w'(s)K aib (t ~ S)ds = Z'(t)g(t) = g'(t). 



(494) 



Now z(t) and g(t) are the restrictions to the real line of entire functions 
of exponential type; so 

g'(r) 



g(r) 



= z'(t). (495) 



Hence g{r) is zero free in the entire plane and is, therefore, of the form 
(Theorem 2.7.1, Ref. 6) 

g(r) = Ae ix \ (496) 

Hence from (495) and (496), z'(t) = i\, and since z (t) is bounded on the 
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real axis it follows that \ = 0; i.e., z(t) = constant and hence 

w(t) = constant. (497) 

APPENDIX E 
Proof of Theorem 6 

We are given 

w(t) = log [z(t)\ (498) 

where the Fourier transform of z(t) vanishes outside [0,Q] and for some 
positive € 

\z(t + iu)\>e forw>0 

-« < t < « (499) 

Thus w(t) is bounded and analytic in the uhp; so by Theorem 1, the 
Fourier transform of w(t) vanishes over (— °°,0). Also we are given that 
the Fourier transform of w(t) vanishes over (a, b) where < a < b 
and 

b - a > fi (500) 

and wish to show that z(t) = constant. 

We proceed as in the proof of Theorem 5 and write 

w(t) = g(t) + h(t) (501) 

where the Fourier transform of g(t) vanishes outside [0,a] and the 
Fourier transform of h(t) vanishes over (— °°,b). We have 

z'(t) 
w'(t)=-f-!- = g'(t) + h'(t) (502) 

or 

z'(t) = z(t)g'(t) + z(t)h'(t). (503) 

Now the Fourier transform of h'(t) vanishes over (— °°,b) so by Corollary 
1 of Theorem 1 the Fourier transform of z(t)h'(t) also vanishes over 
(-oo,fe). By Corollary 2 of Theorem 2, the Fourier transform of z(t)g'(t) 
vanishes outside [0,a + fl]. Since a + Q < b we conclude as in the proof 
of Theorem 5 that 

z'(t)=z(t)g'(t), org'{t) = Z -^- (504) 

and hence that z(t) = Ae ixt and since g{t) is bounded, X = 0. There- 
fore 

z(t) = constant. (505) 
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APPENDIX F 

Lower Bound on the Degree of Certain Polynomials 

Suppose a polynomial of degree v is of the form 

P v (z) = 1 + a n z n + a n+ iz n+1 + ...a v z" (506) 

where \a n \ > 1, and P v {z) is zero free for \z\ < 1. Then 

v > In. (507) 

To prove this assertion we assume 

v < 2n. (508) 

Then assuming that P v (z) is zero free for \z\ < 1, the function 

B, + «,-!« + ...«,«-'' + «' 

1 + o n z n + a n+ i2" +1 + . . . + a^" 
is analytic for | z \ < 1 and 

\f(e i6 )\ =1, -7T < < 7T. (510) 

Then 

f(z)-±b k z k , \z\<\ (511) 



where 

&*=— C* f(e ie )e- ike d6. (512) 

27T %) —TC 

Thus from (510) and (512), 

|M<1. (513) 

However, with the assumption c < 2n we see from (509) that 

b k = a„- k for < k < v - n < n. (514) 

In particular 

b v - n = a n . (515) 

But \a n | > 1, so (515) contradicts (513) and therefore (508) is false. 
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